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Abstract. These lectures explain that comparisons between experiment and theory can 
expose the impact of running couplings and masses on hadron observables and thereby aid 
materially in charting the momentum dependence of the interaction that underlies strong- 
interaction dynamics. The series begins with a primer on continuum QCD, which introduces 
some of the basic ideas necessary in order to understand the use of Schwinger functions as 
a nonperturbative tool in hadron physics. It continues with a discussion of confinement 
and dynamical symmetry breaking (DCSB) in the Standard Model, and the impact of these 
phenomena on our understanding of condensates, the parton structure of hadrons, and the 
pion electromagnetic form factor. The final lecture treats the problem of grand unification; 
namely, the contemporary use of Schwinger functions as a symmetry-preserving tool for the 
unified explanation and prediction of the properties of both mesons and baryons. It reveals that 
DCSB drives the formation of diquark clusters in baryons and sketches a picture of baryons 
as bound-states with Borromean character. Planned experiments are capable of validating the 
perspectives outlined in these lectures. 


1. Continuum QCD: A Primer 

1.1. Quantum chromodynamics 

Quantum chromodynamics (QCD) holds many fascinations, which are somehow hidden in the 
seemingly simple local Lagrangian density that defines the theory: 


-Cqcd = qi[il^[D^]ij - M5ij]qj - 

= qi[il^d^ - m]qi - gG^^q^-f^^^qj - |G“^G^^ 
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( 1 . 2 ) 


where {T“|o = 1,...,8} are the generators of 5C/c(3) in the fundamental representation; G 
represents the gluon gauge-fields; q, the quark matter-fields; g is the coupling constant; and the 
gluon field-strength tensor is 



(1.3) 


A first observation is that QCD works: there is no confirmed breakdown over an enormous 
energy domain: 0 < .E < 8TeV. Furthermore, as these notes will explain, it is probable that 
QCD is a self-contained, nonperturbatively renormalisable and hence well-defined quantum held 





Figure 1. Left panel - Running coupling in QED; and right panel - running coupling in QCD. N.B. 
On the same momentum domain, the QCD coupling changes 500 000-times more than the QED coupling 
and runs in the opposite direction. 


theoryQ If that is so, then QCD is a true theory, not an effective theory. This possibility 
and the lack of any evidence for supersymmetry in experiments at the large hadron collider 
(LHC) have led to a resurgence of interest in the possibility that any reasonable extension of 
the Standard Model will be based on the paradigm established by QCD. One possibility is 
Extended Technicolour, in which electroweak symmetry breaks via a fermion bilinear operator 
in a strongly-interacting non-Abelian theory ii and the Higgs sector of the Standard Model 
becomes an effective description of a more fundamental fermionic theory, similar to the Ginzburg- 
Landau theory of superconductivity. These features and possibilities lend urgency to the problem 
of solving QCD, and hence also to modern programs in hadron physics. 

The underlined term in Eq. (11.31) is the source of everything that is remarkable about QCD 
because it generates gauge-field self-interactions, which have extraordinary consequences. This 
is readily elucidated through a comparison with QED. A characteristic feature of QED is that 
gauge-boson self interactions do not occur at tree level, i.e. at lowest order in perturbation theory. 
The leading contribution occurs at order a^\ and since a ~ 1/137, such interactions take place 
with extremely small probability. In contrast, the non-Abelian character of QCD produces 
tree-level interactions between gauge bosons - the three- and four-gluon vertices generated by 
Gfti/Ga" in Eq. (|l.ip . One might guess that these terms should have a big impact; and doing 
better than guessing led to the 2004 Nobel Prize in Physics [^-[^. 

Relativistic quantum gauge-field theories are typified by the feature that nothing is constant. 
The distribution of charge and mass, the number of particles, etc.; indeed, all the things that 
quantum mechanics holds hxed come to depend upon the wavelength of the tool being used to 
measure them. Couplings and masses are renormalised via processes involving virtual-particles. 
Such effects make these quantities depend on the energy scale at which one observes them, and 
produce the running couplings illustrated in Eig.lH At one-loop order, the QED coupling is 


aQED(Q) =- 7 ^ -py- 

37r Ule 


(1.4) 


where a and nie are renormalised such that they take their observed values when measured at 
the physical electron mass, Q = mg. The negative sign in the denominator signals that fermions 

^ This would make QCD unique. In contrast, it is certainly not true of quantum electrodynamics (QED), 
an Abelian gauge theory, which is perturbatively simple but, alone, undefined nonperturbatively: four-fermion 
operators become relevant in strong-coupling QED and must be included in order to obtain a well-defined 
continuum limit OH. 





screen electric charge. In QCD, the analogous one-loop quantity is 


aQCD(Q) =- Q —, (1-5) 

(33-21V^)ln^ 

AqcD 

where IV/ is the number of quark flavours that are active in the process under consideration and 
^QCD ~ 300 MeV is a mass-scale0 Plainly, quarks screen the colour charge. However, owing 
to their self-interactions, gluons antiscreen this charge. The enormous impact of that effect is 
evident in Fig.[TJ 

The running coupling depicted in the right panel of Fig. [T] can be translated into an interac¬ 
tion between quarks and gluons that depends strongly on separation, i.e. an interaction which 
grows rapidly with separation between the sources. Taking the experimental value of Aqcd; the 
coupling is found to be huge at a separation r ~ where Vp is the proton’s charge radius. 
This is a peculiar circumstance, viz. an interaction that becomes stronger as the participants 
try to separate. It leads one to explore some fascinating possibilities: if the coupling grows so 
strongly with separation, then perhaps it is unbounded; and conceivably it would require an 
inhnite amount of energy in order to extract a quark or gluon from the interior of a hadron? 
This manner of thinking has led to the 

Confinement Hypothesis: Colour-charged particles cannot he isolated and therefore 
cannot be directly observed. They clump together in colour-neutral bound-states. 

This is an empirical fact; but there is no mathematical demonstration. Possibly in the hope 
of focusing attention in that direction, at the turn of the Millennium the Clay Mathematics 
Institute offered a prize of $l-million for a proof that SUfi^) gauge theory is mathematically 
well-dehned (I^ , one necessary consequence of which will be an answer to the question of whether 
or not the confinement conjecture is correct. I will return to this issue. 


13l |. refers to a class of subatomic particles that are 


1.2. Hadron Physics 
The term “hadron,” hrst used in 1962 
composed of quarks and/or gluons and take part in the strong interaction. Well known examples 
are the proton, neutron and pion. The world of hadrons is subdivided into baryons, which are 
hadrons with half-integer spin that contain three valence-quarks; and mesons, hadrons with 
integer spin that contain a valence-quark and -antiquark. Hadron physics is unique at the 
cutting edge of modern science because Nature has provided us with just one instance of a 
fundamental strongly-interacting theory, viz. QCD: basic science has never before confronted 
such a challenge as solving this theory. 

In facing this task, hadron physics has developed into an international research endeavour 
of remarkable scope. Indeed, before the end of this decade the field will be operating a host of 
upgraded or new facilities and detectors. An illustrative list may readily be compiled: Beijing’s 
electron-positron collider; J-PARC, the Japan Proton Accelerator Research Complex, 150km 
NE of Tokyo; the ALICE and COMPASS detectors at CERN; and, in the USA, both RHIC 
(Relativistic Heavy Ion Collider), at Brookhaven National Laboratory, and JLab (Thomas 
Jefferson National Accelerator Eacility), in Newport News, VA. RHIC focuses primarily on the 
strong-interaction phase transition, i.e. physics just 10/is after the Big Bang; and JLab explores 
the nature of cold hadronic matter. Regarding the upgrade of JLab [1^, commissioning of Hall- 
A and Hall-D are now taking place, a process that is itself expected to yield new physics results, 
and project completion is expected in 2017. In addition to these existing facilities and those 


^ Aqcd is the natural mass-scale of QCD. It is dynamically general^ via quantisation. The appearance of Aqcd 
spoils the conformal invariance of the classical massless theory [9l4ll||. The value of Aqcd must be determined 
empirically within the Standard Model; and it sets the scale for all dynamically generated masses described herein. 





under construction, excitement is also generated by the discovery potential of proton-nucleus 
collisions at the LHC [I^ and new machines, such as an electron ion collider (EIC) (l^ . 

A theoretical understanding of the phenomena of hadron physics requires use of the full 
machi nery of relativistic quantum field theory. Based on the relativistic dynamics elucidated by 
Dirac |17l ]. it is the only known way to reconcile quantum mechanics with special relativity. 

It is worth highlighting that the unification of relativity with quantum mechanics took some 
time; and questions still remain as to a practical implementation of a Hamiltonian formulation 
of the relativistic quantum mechanics of interacting systems. The Poincare group has ten 
generators: six associated with the Lorentz transformations (rotations and boosts); and four 
associated with translations. Quantum mechanics describes the time evolution of a system with 
interactions and that evolution is generated by the Hamiltonian, or some generalisation thereof. 
However, the Hamiltonian is one of the generators of the Poincare group, and it is apparent 
from the Poincare algebra that boosts do not generally commute with the Hamiltonian. This is 
partly because the existence of antiparticles is concomitant with relativistic quantum mechanics, 
i.e. the equations of relativistic quantum mechanics admit negative energy solutions. However, 
once one allows for negative energy, then particle number conservation is lost: 


Egystem — Egyatejii -\- {^Ep^ -\- Ep^'j T * * * ad infinitum^ 


( 1 . 6 ) 


where Ep = —Ep. This poses a fundamental problem for relativistic quantum mechanics: few 
particle systems can be studied, but the study of (infinitely) many bodies is difficult and no 
general theory currently exists. Following from these observations, the state vector calculated 
in one momentum frame will not be kinematically related to the state vector for the original 
system in another frame, a fact that makes a new calculation necessary in every momentum 
frame. The discussion of scattering, which takes a state of momentum p to a different state with 
momentum p' is therefore problematic [iS-Q. 

Relativistic quantum field theory provides a way forward. In this framework the fundamental 
entities are fields, which can simultaneously represent infinitely many particles. The neutral 
scalar held, (p{x), provides an example. One may write 



1 

(27r)3 2ujk 


a(A:)e-*^'^ + ot(A:)e*^-^ 


(1.7) 


where: Uk = y 1^1^ + is the relativistic dispersion relation for a massive particle; the four- 

vector (k^) = {ujk,k)] a{k) is an annihilation (creation) operator for a particle (antiparticle) 
with four-momentum k {—k)] and (k) is a creation (annihilation) operator for a particle 
(antiparticle) with four-momentum k {—k). With this plane-wave expansion of the held one 
may proceed to develop a framework in which the nonconservation of particle number is not a 
problem. That is crucial because key observable phenomena in hadron physics are essentially 
connected with the existence of virtual particles. 

Relativistic quantum held theory has its own problems, however. For example, the question 
of whether a given quantum held theory is rigorously well dehned is an unsolved mathematical 
problem. All relativistic quantum held theories admit analysis via perturbation theory, and 
perturbative renormalisation is a well-dehned procedure that has long been used in QFD and 
QCD. However, the rigorous dehnition of a theory means proving that the theory makes sense 
nonperturbatively. This is equivalent to proving that all the theory’s renormalisation constants 
are nonperturbatively well-behaved. 

Understanding the properties of hadrons requires solving QCD. However, despite the 
possibility highlighted above and discussed further below, QCD is not now known to be a 
rigorously well-dehned theory. It might not have a solution. Thus, more properly, the task of 
hadron physics is to determine whether QCD is truly the theory of the strong interaction. 






Physics is an experimental science and experiment advances quickly. Hence, developing an 
understanding of observable phenomena cannot wait on mathematical rigour. Assumptions must 
be made and their consequences explored. Hadron physics practitioners therefore assume that 
QCD is (somehow) well-defined and follow where that idea may lead. This means exploring 
and mapping the hadron physics landscape with well-understood probes, such as the electron at 
JLab; and employing established mathematical tools, refining and inventing others, in order to 
use the Lagrangian density of QCD to predict what should be observable real-world phenomena. 

A primary aim of the world’s current hadron physics programmes in experiment and theory 
is to determine whether there are any contradictions with what can actually be proved in QCD. 
There are none today; but that doesnt mean there are neither puzzles nor controversies. The 
interplay between experiment and theory is the engine of discovery and progress: the discovery 
potential of both is high. Much has been learnt in the last ten years; and these lectures will 
describe some of those discoveries and provide a perspective on their meaning. Many of the 
most important questions in basic science are the purview of Hadron Physics. 


1.3. Green Funetions and Propagators 

The Dirac equation provides the starting point for a Lagrangian formulation of the quantum 
field theory for fermions interacting via gauge boson exchange. For a free fermion with mass m, 

[i(^x -m]'ijj{x) = 0, (1.8) 

where ijj^x) is the fermion’s “spinor” - a four component column vector, while in the presence 
of an external electromagnetic field the fermion’s wave function obeys 


[i(^x — eJ^ix) — m] iplx) = 0 . 


(1.9) 


which is obtained, as usual, via “minimal substitution:” ^ — eA^ in Eq. (I1.8jl . These 

equations have a manifestly covariant appearance and a proof of covariance is given in the early 
chapters of Refs. 21. 1^. 

The Dirac equation is a partial differential equation. A general method for solving such 
equations is to use a Green function, which is the inverse of the differential operator that 
appears in the equation. The analogy with matrix equations is obvious and can be exploited 
heuristically. 

Equation (11.91) yields the wave function for a fermion in an external electromagnetic field. In 
this connection, consider the operator obtained as a solution of the following equation 


— eJjL{x') — m] S{x',x) = 1 S‘^{x' — x). 


( 1 . 10 ) 


It is immediately apparent that if, at a given spacetime point x, ip{x) is a solution of Eq. (11.91) . 
then 

ijj(x^) := j d^x S{x\x)'ij}{x) (1-11) 

is a solution of 

[i(j)x' — eJ^ix') — m]'i/'(x') = 0 , (1-12) 

i.e. S{x',x) has propagated the solution at x to the point x'. 

This effect is equivalent to the application of Huygens’ principle in wave mechanics: if the 
wave function at x, 'il’{x), is known, then the wave function at x' is obtained by considering ip(x) 
as a source of spherical waves that propagate outward from x. The amplitude of the wave at x' 
is proportional to the amplitude of the original wave, the constant of proportionality 

is the propagator (Green function), S{x',x). The total amplitude of the wave at x' is the sum 
over all points on the wavefront, viz. Eq. (II.lip . 





Such an approach is practical because all physically reasonable external fields can only be 
nonzero on a compact subdomain of spacetime. Therefore the solution of the complete equation 
is transformed into solving for the Green function, which can then be used to propagate the free- 
particle solution, already found, to arbitrary spacetime points. However, obtaining the exact 
form of S{x',x) is impossible for all but the simplest cases (see, e.g. Ref. [H, [^). 

In the absence of an external field Eq. (jl.lOl) becomes 

— m] S{x',x) = 1 d^{x'— x). (1-13) 

Assume a solution of the form: 

Soix',x) = Soix'-x) = I (1.14) 

so that substituting yields 

- m) So(p) = 1; z.e. So(p) = ^ ^2 ■ 

To obtain the result in configuration space one must adopt a prescription for handling the 
on-shell singularities in S{p), viz. the pole at = m?, and that convention is linked with the 
boundary conditions applied to Eq. (jl.ldp . An obvious and physically sensible definition of the 
Green function is that it should propagate positive-energy-fermions and -antifermions forward 
in time but not backwards in time, and vice versa for negative energy states. 

Gonsider now that the wave function for a positive energy free-fermion is 

'iIj^^\x) = ti(p) . (1-16) 

The wave function for a positive-energy antifermion is the charge-conjugate of the negative- 
energy fermion solution: 


^(+)(x) = C7° , (1.17) 

where C = and (•)"’" denotes matrix transposed It is thus evident that the physically 

sensible Sq{x' — x) must only contain positive-frequency components for x'q — xq > 0, i.e. it must 
be proportional to the positive-energy projection operator A+(p) = {i) + m)/(2m). One can 
ensure this via a small modification of the denominator in Eq. (I1.15p : 


So{p) 


if) + m 

p2 _ jyy2 


^ + m 

p2 _ yyy2 _|_ ’ 


(1.18) 


with ry —)• 0'*' at the end of all calculations. Inserting Eq. (11.1811 into Eq. (I1.14p is equivalent 
to evaluating the p^ integral by employing a contour in the complex-p*^ plane that is below 
the real-p® axis for p^ < 0, and above it for p^ > 0. This prescription defines the Feynman 
propagator. 

Let’s now return to Eq. (II.10^ : 


— eJ^{x'^ — m] S{x',x) = 1 5^{x'— x), (1-19) 

® This defines the operation of charge conjugation. N.B. The form of C introduced here is appropriate to a 
particular Dirac matrix representation. Nevertheless, the unitary equivalence of Clifford algebra representations 
ensures that generality is not lost [iH] . 


















which dehnes the Green function for a fermion in an external electromagnetic field. As 
mentioned, a closed form solution of this equation is impossible in all but the simplest held 
conhgurations. Is there, nevertheless, a way to construct an approximate solution that can be 
systematically improved? 

To achieve that one rewrites the equation thus: 

—'m\S{x\x) = — x) ^ Six'^x) ^ (1.20) 

which, as easily seen by substitution, is solved by 

Six ,x) = Soix'- x) + e j d'^ySoix'-y)4iy)Siy,x) (1.21) 

= Soix' -x) + e J d‘^y 5o(x' - y)4iy) Soiy - x) 

+e^ f d'^yi [ d'^y2SQix' -yi)4.iyi)Soiyi-y2)4{y2)Soiy2-x) 

+ [...]. ( 1 . 22 ) 

This perturbative expansion of the full propagator in terms of the free propagator provides 
an archetype for perturbation theory in quantum held theory. One obvious application is the 
scattering of an electron/positron by a Coulomb held, which is an example explored in Sec. 2.5.3 
of Ref. [2^. Equation (|1.22l) is a hrst example of a Dyson-Schwinger equation (DSE) (^ . 

In quantum held theory this Green function has the following interpretation: 

(i) It creates a positive energy fermion (antifermion) at spacetime point x; 

(ii) Propagates the fermion to spacetime point x', i.e. forward in time; 

(iii) Annihilates this fermion at x'. 

The process can equally well be viewed as 

(i) The creation of a negative energy antifermion (fermion) at spacetime point x'] 

(ii) Propagation of the antifermion to the spacetime point x, i.e. backward in time; 

(iii) Annihilation of this antifermion at x. 

Other propagators have similar interpretations. 


l.J^. Dyson Schwinger Equations: e.g. Photon Vacuum Polarisation 

Local gauge theories are the keystone of contemporary hadron and high-energy physics. Such 
theories are difficult to quantise because the gauge dependence is an extra non-dynamical degree 
of freedom that must be handled. The modern approach is to quantise the theories using the 
method of functional integrals, attributed to Eeynman and Kac. References 


11 


provide 

useful overviews of the technique, which replaces canonical second-quantisation. 

Beginning with the field equations of quantum field theory, it has long been known that one 
can derive a system of coupled integral equations interrelating all of a theory’s Green functions 
27]. This collection of a countable infinity of equations is called the complex of DSEs. It is an 


intrinsically nonperturbative complex, which is vitally important in proving the renormalisability 
of quantum field theories, and, at its simplest level, the complex provides a generating tool for 
perturbation theory. I will illustrate a nonperturbative derivation of one equation in this complex 
within the context of quantum electrodynamics QED. The derivation of others follows the same 
pattern. 







1.4-1- Generating Functional and Maxwell’s Equations. Let’s begin with the action for QED 
with Nf flavours of electromagnetically active fermions: 




/ 


d‘^x 


Nf 

/=i 






fiiy 


_ 


2Ao 


d^Af,{x) d‘'Ay{x) 


(1.23) 

In this manifestly Poincare covariant action: [x), {x) are the elements of the Grassmann 

algebra that describe the fermion degrees of freedom, so that, e.g. 'df,g, = 0; ml 

are the fermions’ bare masses and el their bare charges; and j4^(x) describes the gauge boson 
[photon] field, with 


= df,Ay - dyAf,, (1.24) 

and Ao is the bare Lorentz gauge fixing parameter. [N.B. To describe an electron, the physical 
charge Cf < 0.) 

The generating QED functional is defined via the action in Eq. (I1.23P : 


= JpA^] expjz I 


d^x 


- F,,{x) - ^ d^A,{x) d^'A^ix) 


+ ^ -ml + el^l^ 4)^ + J^{x)Afj,{x) + {x)'ijj^ (x) + {x)^^ {x) \ , 

/=! -I J 

(1.25) 


where is an external source for the electromagnetic field, and are external sources for 

the fermion fields that are also elements in the Grassmann algebra. An n*'^-order functional 
derivative of with respect to any combination of its arguments yields an expectation 

value of a product of fields, i.e. a moment of the measure. It follows that so long as the integral 
is well defined, then all the theory’s physical content is expressed in Eq. (|1.25p : but is the integral 
well defined? 

The product [P'0] in Eq. ()1.25l) indicates that one is supposed to integrate over 

all possible values of all the fields at all the spacetime point in the Universe. A priori, that 
would seem an impossible task. Gonsequently, any meaningful definition of relativistic quantum 
field theory is formulated in Euclidean space. The oscillating weight factor then becomes an 
exponential damping factor, which acts to suppress contributions from all configurations that are 
distant from the action’s stationary point(s). Eor the remainder of this section I will nevertheless 
continue to work in Minkowski space, ignoring the problem because most readers will be more 
familiar with Minkowski space conventions. Nothing is lost in doing this: all the manipulations 
performed in this sketch of the derivation of the DSE for the photon vacuum polarisation will 
be purely formal, with no attempt at mathematical rigour. 

The derivation of a DSE follows simply from the observation that the integral of a total 
derivative vanishes, viz. 

/ 

given appropriate boundary conditions. Gonsider, therefore, 
r _ (5 i(S[Au,'4’,'4’]+f 


0 
















X exp 


i ( + f <fx 


5S 


5A^{x) 


5 _ 5_ 

iSJ' iS^' i6^ 


+ Mx)}w[j^,^,(], 


(1.27) 


where the last line has meaning as a functional differential operator acting on the generating 
functional. It may be interpreted as an implementation in QFT of the “principle of stationary 
action”, where the generating functional replaces the classical action, i.e. as an Euler-Lagrange 
equation in QFT. 

To proceed from here it is advantageous to introduce the generating functional of connected 
Green functions, viz. defined via 


Then, making use of the fact that differentiating Eq. (|1.23l) gives 


6S 


dA^{x) 

Eq. (|1.27p becomes 


dpdPg^y - ( 1 - ^ ) dpdy 




+ ^^4 '>p4x)'ypipf {x ), 
/ 


6Z 


(1.28) 


(1.29) 


^ f f 5Z 5Z 5 \ 6iZ 


where I have divided through by Equation (ll.3Up represents a compact form of the 

nonperturbative equivalent of Maxwell’s equations. 


1 . 4 . 2 . One Particle Irreducible Green Functions. The next step is to introduce the generating 
functional for one-particle-irreducible (IPI) Green functions: T)^^, ij}, ip], which is obtained from 
Z[J^,^,^] via a Legendre transformation 


z[jp,c,c] = n\,i’,4 + 


I 


dP^x 




(1.31) 


A one-particle-irreducible n-point function or “proper vertex” contains no contributions that 
become disconnected when a single connected m-point Green function is removed, e.g. via 
functional differentiation. This is equivalent to the statement that no diagram representing 
or contributing to a given proper vertex separates into two disconnected diagrams if only one 
connected propagator is cut. (A detailed explanation is provided in Ref. [2^, pp. 289-294.) 

The following identities are plain from the definition of the generating functional, Eq. (11.251) : 


6Z 

6J>^{x) 


\{x) , 


6Z 


^p{x) , 


6Z 


-Tpix) , 


(1.32) 


where here the external sources are nonzero. Hence T in Eq. (I1.31D must satisfy 


6T 

6Ai^{x) 




.5r 

d-ipl (x) 


-4ix) , 


5T 

S-ipf (x) 


= 4ix) ■ 


(1.33) 





































Notably, since the sources are not zero: 

XA (^\ 

A,{x)=A,ix-,[J„C,C]) => ( 1 - 34 ) 

with analogous statements for the Grassmannian functional derivatives. 

It is easy to see that setting ip = Q = ijj after differentiating T gives zero unless there are 
equal numbers of ■0 and 0 derivatives: any integral that is odd under 0 —)■ —0, 0 —0 must 

vanish. Similarly, with = 0 the generating functional for IPI Green functions is even under 
Ajj^ —>• —A^, from which follows Furry’s Theorem, i.e. the vacuum expectation value of any odd 
number of electromagnetic currents is zero. 

Gonsider the operator and matrix product (with spinor labels r, s, t) 
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- (Tz- 


6‘^Z 




^^r(xXHz) <50n2:)00y) 


0 


0=0 
0 = 0 


which, using Eqs. (I1.32jl . (I1.33jl . simplifies as follows: 


[a 

J 6 Cf{x)S^PH^) 


0 


0 = 0 
0 = 0 


(y) 



0 = 0 


= 5rs - y). 


(1.35) 


(1.36) 


Returning to Eq. (I1.30p and setting C = 0 = ? one obtains 


sr 


6At^{x) 


ll;='lp=0 


dpdPguu - ( 1 - ) dad> 


after making the identification 
S^{x,y; [Afj]) = - 


Ao ' " 

5^Z 


-4''(a:) - i ^e^tr 7^S'^(a:,x; [A/,]) , (1.37) 
/ 


5^Z 


^CHyWix) SCHx)C^{y) 


(no summation on /), 


(1.38) 


which is the connected Green function that describes the propagation of a fermion with flavour 
/ in an external electromagnetic field (Recall the free fermion Green function in Eq. (Il.ldp .i 
Notably, as a direct consequence of Eq. (I1.35P the inverse of this Green function is given by 


S^{x,y;[A]) ^ = 


<52 r 


{x)5'ip'f {y) 


(1.39) 


xlj='ip=0 


It is a general property that such functional derivatives of the generating functional for IPI 
Green functions are related to the associated propagator’s inverse. Glearly the vacuum fermion 
propagator or connected fermion 2-point function is 


S^{x,y) := S^{x,y; [A^ = 0]). 

Such vacuum Green functions are keystones in quantum field theory. 

At this point one differentiates Eq. (| 1.3711 with respect to Ay{y) and sets Jp{x) = 0: 

d^p 


(1.40) 


6Ai^{x)6A'^{y) 


Ap = 0 
-0 = 0 = 0 


dpdPgay - 1 - 
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Ao 


dpdy 


S {x - y) 


-i^e^tr 
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In 


52 p 


5Av{y) \ 5'4)f {x)5il)f {x) 


-1 


1^='0 = O i 


(1.41) 










































The left-hand-side (Ihs) is readily understood. Just as Eqs. (I1.39jl . (11.4011 define the inverse 
of the fermion propagator, here one has 


{D-^r{x,y) := 




dAt^{x)dA'^{y) 


A^ = 0 

Ip = Ip = 0 


(1.42) 


viz. the inverse of the photon propagator. 

On the other hand, the rhs must be simplihed and interpreted. First observe that 




6Ay{y) \ 6ipf{x)6'ipf(x 


-1 


= — J d!^ud^w 


^=•0=0 y 


5ipf {x)6'ipf {w) 


-1 






'0='0=O y 


5Ay{y) 5ipf {u)6'ipf {w) I Sipf {w)6ip-f (x) 


= O / 
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which is an analogue of the result for finite dimensional matrices: 

Equation (ll.43p involves the IPI 3-point function (no summation on /) 

This is the proper fermion-gauge-boson vertex. At leading order in perturbation theory 

r^(x, y; z) = -fu 5‘^{x - z) <5^(y - z ), (1.46) 

a result which can be obtained via explicit calculation of the functional derivatives in Eq. (I1.45p . 

1 . 4 . 3 . Photon Self Energy. Defining the gauge-boson vacuum polarisation: 

Ilfj,^{x,y) = i'^ielf f d!^zid!^Z2\.T {x,zi)Tl{zi,Z2]y)S^{z2,x) , (1.47) 

/ 

which is illustrated in Fig.[2l it is immediately apparent that Eq. (|1.4ip may be expressed as 
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(1.48) 


In general, the gauge-boson vacuum polarisation, or “photon self-energy,” describes the 
modification of the gauge-boson’s propagation characteristics owing to the presence of virtual 
particle-antiparticle pairs in quantum field theory. In particular, the photon vacuum polarisation 
is an important element in the description of many physical processes, for instance —?> e'’'e“. 

The propagator for a free gauge boson is obtained when = 0 and may be written, in 
momentum space: 
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Figure 2. Upper panel ~ photon vacuum 
polarisation in Eq. (11.471) . Lower panel - DSE for 
the photon propagator, whose kernel is the photon 
vacuum polarisation. 


In the presence of interactions, i.e. for 11^,^ 7 ^ 0 in Eq. (I1.48jl . this becomes 

—g^'^ + {q^q''/[g^ + ig]) i ^ 9^'?" 
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where I have used the Ward-Green-Takahashi identity [2^-[4( 
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which means that one can write 


n^"(<?) = + <?V) n(g2). (1.52) 

n(q^) is the polarisation scalar. In QED, it is independent of the gauge parameter, Aq- The 
choice Ao = 1 is called “Feynman gauge”; and it is useful in perturbative calculations because it 
simplifies the n(g^) = 0 gauge boson propagator enormously. In nonperturbative applications, 
however, Aq = 0 , “Landau gauge,” is most useful because, inter alia, it ensures that the gauge 
boson propagator is itself transverse. Equation (II.50^ corresponds to the DSE expressed in the 
lower panel of Fig. [21 

The longitudinal Ward-Green-Takahashi identities (WGTIs) (28l-[^ are relations satisfied by 
n-point Green functions, relations which are an essential consequence of a theory’s local gauge 
invariance, i.e. local current conservation. They can be proved directly from the generating 
functional and have physical implications. For example, Eq. (|1.52p ensures that the photon 
remains massless in the presence of charged fermions @ 

In the absence of external sources for fermions and gauge bosons, Eq. (|1.47p can easily be 
represented in momentum space, for then the 2-and 3-point functions appearing therein must 
be translationally invariant and hence they can be expressed in terms of Fourier amplitudes, i.e. 
one has 


in^^(g) = - j -^^ii[{ig^,){iS^{l)){iTf + q)){iS{l + q))]. (1.53) 


It is the reduction to a single integral that makes momentum space representations most widely 
used in continuum calculations. 

The vacuum polarisation in QED is directly related to the running coupling constant, which 
is a connection that makes its importance obvious. The connection is not so direct in QGD; but, 
nevertheless, the polarisation scalar is a key component in the evaluation of the strong running 
coupling. 

* A discussion of longitudinal WGTIs can be found in Ref. [^, pp. 407-411, and Ref. [dH , p p. 299-303; and their 
generalisation to non-Abelian theories as “Slavnov-Taylor” identities is described in Ref. ]4|, Chap. 2. There are 
also transverse WGTIs [3ll-l4o|. which have additional, important consequences. 
















In the above analysis we saw that second derivatives of the generating functional, r[^^, -0, tp], 
give the inverse-fermion and -photon propagators and that the third derivative gave the proper 
photon-fermion vertex. In general, all derivatives of this generating functional, higher than 
two, produce a corresponding proper vertex, where the number and type of derivatives give the 
number and type of proper Green functions that it can serve to connect. 

2. The Ins and Outs of Mesons 

2.1. Elegance Is Not Everything 

It is worth reiterating at this point that physics is an empirical science. I stress this because, faced 
with challenges in applying fundamental theories to the observed Universe, some researchers have 
called for a change in how theoretical physics is done [i^. They have begun to argue that if 
a theory is sufficiently elegant, it need not be tested experimentally. Prominent amongst the 
“elegance will suffice” advocates are some string theorists and cosmologists. This is plainly 
because their elegant theories have not found experimental confirmation at the LHC. This is 
not science; and it brings me to recollect the words of perhaps the world’s most famous PI (i^ : 

I have no data yet. It is a capital mistake to theorize before one has data. Insensibly one 
begins to twist facts to suit theories, instead of theories to suit facts. 

Thus, once more, physics is an empirical science. If a claim cannot be tested empirically, 
then it does no lie within the realm of physics; and no claim can be considered proven unless 
it’s verified experimentally. Viewed from this perspective, the theoretical physicist’s task is to 
develop theories and elicit their observable content, viz. elucidate a diverse array of the theory’s 
most basic measurable predictions. In due time experimental outcomes will then decide which 
facets of Nature are explained by the theory and where lie the boundaries of its applicability. 
Perhaps there is no theory of everything] or, perhaps, the theory of everything will derive its 
elegance by following the paradigm established by QCD? 


2.2. Jefferson Lab 

A context for the remainder of these lectures is provided by the experimental programmes 
underway and planned at JLab in Newport News, Virginia. Funds for the development of 
research plans and designs for the Continuous Electron Beam Accelerator Facility (CEBAF) 
at JLab were initially provided in 1984, and construction began in 1987. Seven years later, in 
1994, the facility achieved its design capability, delivering 4 GeV electron beams on targets in 
the three associated experimental halls. The goal was to “write the book” about the strongest 
known force in Nature - the force that holds nuclei together - and determine how that force can 
be explained in terms of the gluons and quarks of QCD. 

An aim of the original JLab programme was verification of the following prediction: at 
energy-scales greater than some minimum value, Qq, whose value was not determined by the 
proof, hadron exclusive form should behave as follows: 
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where the integer N = Uquarks — I, with Uquarks counting the lowest possible number of quarks 
and/or antiquarks in the target. The power-law term gives rise to parton-model scaling and can 
be explained on dimensional grounds. The distinctive signature of QCD lies in the logarithmic 
factor, the exponent on which, 7 ^, can be computed and whose appearance signals scaling 
violation, which is well documented in deep inelastic scattering [^. It was optimistically 
imagined that Qq might be as small as 1 GeV and so CEBAF was designed to reach 4 GeV. 








In the ten years following achievement of design capacity, numerous fascinating experimental 
results were obtained at JLab, including an empirical demonstration that the distribution of 
charge and magnetization within the proton are completely different 
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where nip is the proton mass. This fact, which overturned a longstanding particle physics 
paradigm, along with a range of related observations, is explained by the presence of strong, 
nonpointlike, electromagnetically-active scalar and axial-vector diquark correlations within the 
nucleon (see Refs. [57-65] and Sec. [3] herein). However, notwithstanding the breadth of JLab’s 
programme and the associated excitement, no experiment has produced an unambiguous signal 
for the realisation of Eq. (12.ip . i.e. as yet we have no clear sign of parton model scaling and 
certainly no sign of scaling violations. Partly owing to this but also with a vast array of new 
experimental tests of the Standard Model’s strong interaction sector in mind 1^, in 2004 it was 
decided that CEBAF should be upgraded so that it could deliver 12 GeV electrons on targets. 
That upgrade is now complete and commissioning beams are being sent to the experimental 
halls. 

The JLab 12 facility and, indeed, an array of modern accelerators worldwide, will confront 
a diverse range of scientihc challenges. In the foreseeable future, we will know the results of a 
search for hybrid and exotic hadrons, the discovery of which would force a dramatic reassessment 
of the distinction between matter and force helds in Nature@ Opportunities provided by new 
data on hadron elastic and transition form factors will be exploited, yielding insights into the 
infrared running of QCD’s coupling and dressed-masses, revealing those correlations which are 
key to hadron structure, exposing the facts or fallacies in contemporary descriptions, and seeking 
verification of Eq. (|2.1I) in numerous processes. Precise experimental studies of the valence- 
quark region will proceed, with the results being used to confront predictions from theoretical 
computations of distribution functions and distribution amplitudes - computation is critical 
here because without it no amount of data will reveal anything about the theory underlying the 
phenomena of strong interaction physics. In addition, the international community will seek and 
exploit opportunities to use precision-QCD as a probe for physics beyond the Standard Model. 
Each of the pieces in this body of exploration, however, can be viewed as steps directed toward 
understanding a single overarching puzzle within the Standard Model, viz. what is confinement 
and how is it related to dynamical chiral symmetry breaking (DCSB) - the origin of the vast 
bulk of visible mass in the Universe? 


2.3. Dyson-Schwinger Equations for Hadron Physics 

In order to match the rate at which progress is being made experimentally and, better, to guide 
and enhance those programmes, flexible, responsive theoretical approaches are needed: methods 
that are capable both of rapidly providing an intuitive understanding of complex problems 
and illuminating a path toward answers and new discoveries. In this milieu, notwithstanding 
its steady progress toward results with input parameters that approximate the real world, the 
numerical simulation of lattice-regularised QCD (IQCD) will not suffice. Approaches formulated 
in the continuum and inspired by, based upon, or connected directly with QCD are necessary. 
Prominent amongst such tools are QCD Sum Rules |^,[^ and DSEs 

These lectures focus on the use of DSEs in cold, sparse hadron physics. (Applications to hot, 
dense QCD may be traced from Ref. [^ .) The DSEs are well suited to the study of QCD because, 

® Exotic mesons are states which possess quantum numbers that are not possible in two-body quark-antiquark 
systems that can be described by quantum mechanics whereas hybrid mesons have quark model quantum numbers 
but unusual decay patterns. Both types of systems are possible in QCD if one admits the concept of “constituent” 
or valence gluon content. 






as remarked above, their simplest application is as a generating tool for perturbation theory. 
Since QCD is asymptotically free, that materially reduces model dependence in intelligent 
applications because the interaction kernel in each DSE is known for all momenta within the 
perturbative domain: k'^ > 2 GeV^. Any model that needs building need then only make 
statements about the long-range behaviour of the kernels. This is good because DSE solutions 
are Schwinger functions, i.e. propagators and vertices; and since all cross-sections are built 
from Schwinger functions, the approach connects observables with the long-range behaviour 
of the theory’s running coupling and masses. Consequently, feedback between theoretical 
predictions and experimental tests can then rehne the statements and lead to an understanding of 
these fundamental quantities. Those predictions are wide-ranging because the DSEs provide a 
nonperturbative, continuum approach to hadron physics and can therefore address questions 
pertaining to, e.g.: the gluon- and quark-structure of hadrons; and the roles of emergent 
phenomena - conhnement and DCSB - and the connections between them. Indeed, the past 
decade has seen the emergence of a novel understanding of gluon and quark confinement and its 
consequences; and we are arriving at a clear picture of how hadron masses emerge dynamically 
in a universe with light quarks, viz. DCSB. Furthermore, computations of ground-state hadron 
wave functions with a direct connection to QCD are now available. They reveal that quark- 
quark correlations are crucial in baryon structure j^; and experimental evidence in support of 
this prediction is accumulating (^ . 


2.4- Confinement 

This notion has already been mentioned. However, in order to consider the concept further it is 
actually crucial to define the subject. That problem is canvassed in Sec. 2.2 of Ref. (t^ ]. which 
explains that the potential between infinitely-h eavy quarks measured in numerical simulations of 
quenched IQCD - the so-called static potential [7^ - is irrelevant to the question of conhnement 
in our Universe, in which light quarks are ubiquitous and the pion is unnaturally light. This is 
because light-particle creation and annihilation effects are essentially nonperturbative in QCD 
and so it is impossible in principle to compute a quantum mechanical potential between two 
light quarks [zi-lTi. This means there is no hux tube in a Universe with light quarks and 
consequently that the hux tube is not the correct paradigm for conhnement. 

DCSB is critical to this property of the Standard Model because it ensures the existence 
of nearly-massless pseudo-Goldstone modes (pions), each constituted from a valence-quark and 
-antiquark whose individual Lagrangian current-quark masses are < 1% of the proton mass (^ . 
In the presence of these modes, no hux tube between a static colour source and sink can have 
a measurable existence. To verify this statement, consider such a tube being stretched between 
a source and sink. The potential energy accumulated within the tube may increase only until 
it reaches that required to produce a particle-antiparticle pair of the theory’s pseudo-Goldstone 
modes. Simulations of IQCD show [ 


77l | that the hux tube then disappears instantaneously 


along its entire length, leaving two isolated colour-singlet systems. The length-scale associated 
with this effect in QCD is ~ (1/3) fm and hence if any such string forms, it would dissolve 
well within a hadron’s interior. 

An alternative perspective associates conhnement with dramatic, dynamically-driven changes 
in the analytic structure of QCD’s propagators and vertices. In this realisation, conhnement is a 
dynamical process. In fact, as will subsequently be explained, contemporary theory predicts that 
both gluons and quarks acquire running mass distributions in QCD, which are large at infrared 
momenta (see, e.g. Refs. [^, The generation of these masses leads to the emergence 

of a length-scale ? ~ 0.5 fm, whose existence and magnitude is evident in all existing studies 
of dressed-gluon and -quark propagators and which characterizes a dramatic change in their 
analytic structure. In models based on such features [1^, once a gluon or quark is produced, it 
begins to propagate in spacetime; but after each “step” of length <j, on average, an interaction 
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Figure 3. Dressed-quark self energy, which is 
the dynamical content of QCD’s most basic fermion 
gap equation. The kernel is composed from the 
dressed-gluon propagator (spring with dark circle) 
and the dressed-quark-gluon vertex (light-circle). 
The equation is nonlinear owing to the appearance 
of the dressed-quark propagator (line with dark 
circle). This image encodes every valid Feynman 
diagram relevant to the quark dressing process 
and, additionally, e.g. any topologically nontrivial 
contributions that might be important. (Momentum 
flows from right-to-left.) 


occurs so that the parton loses its identity, sharing it with others. Finally a cloud of partons is 
produced, which coalesces into colour-singlet final states. Such pictures of parton propagation, 
hadronisation and conhnement can be tested in experiments at modern and planned facilities. 

2.5. Dynamical chiral symmetry breaking 

Whilst the nature and realisation of confinement in empirical QCD is still being explored, DCSB; 

namely, the generation of mass from nothing, is a theoretically-established nonperturbative 

feature of QCD. It is important to insist on the term “dynamical,” as distinct from spontaneous, 

because nothing is added to QCD in order to effect this remarkable outcome and there is no 

simple change of variables in the QCD action that will make it apparent. Instead, through the 

act of quantising the classical chromodynamics of massless gluons and quarks, a large mass-scale 

is generated. DCSB is the most important mass generating mechanism for visible matter in the 

Universe, being responsible for approximately 98% of the proton’s massU 

A fundamental expression of DCSB is the behaviour of the quark mass-function, M(p), which 

fyl ! \± / < 

is a basic element in the dressed-quark propagatoiU 

S{p) = l/[i-^ ■ pA{p‘^) + B{j?)] = Z{p‘^)/[i-fp + M(jl ^)], (2.3) 

which may be obtained as a solution to QCD’s most basic fermion gap equation (see Fig.[3]). The 
nontrivial behaviour of the mass function, depicted in Fig. 01 arises primarily because a dense 
cloud of gluons comes to clothe a low-momentum quark; and explains how an almost-massless 
parton-like quark at high energies transforms, at low energies, into a constituent-like quark, 
which possesses an effective “spectrum mass” Mq ~ 350 MeV. Consequently, the proton’s mass 
is two orders-of-magnitude larger than the sum of the current-masses of its three valence-quarks. 

One might ask just how the self-energy depicted in Fig.[3] is capable of generating mass 
from nothing, viz. the m = 0 curve in Fig. 01 which cannot arise in the classical theory. 
The answer lies in the fact that Fig.O is a decratively simply picture. At the very least, it 
corresponds to a countable infinity of diagramsjj all of which can potentially contribute. To 

® This means that the Higgs boson plays almost no role in generating the mass of anyone reading this text. 
Here “almost” means that there must be something like a Higgs mechanism in the Standard Model because at 
least one of the light quarks should have a nonzero mass so that the pion, too, acquires a nonzero mass, with 
magnitude < Aqcd. In the presence of an infinite-range strongly-interacting boson, the Universe could not exist 
in its present form. 

^ I have now switched to a Euclidean metric, in which the Dirac matrices are defined via {7^, 7 i/} = 25^,/, 7), = 7^, 
75 = 74717273 , tr[757^7,.7p7a] = cTp,. = (i/ 2 )[ 7 p, 7 ,/]; and a-b = Y.t=i Pp timelike =7 < 0. 

* In principle, the dressed-quark-gluon vertex involves terms that express DCSB and consequently contains 
information about any topologically nontrivial configurations that are important in this connection |8(1 . 
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Figure 4. Dressed-quark mass function, M{p) 
in Eq. (12.31) : solid curves ~ DSE results, explained 
in Refs. jSflI. 1^ . “data” - numerical simulations of 
lattice-regularised QCD [s^- {N.B. m = 70MeV 

is the uppermost curve and current-quark mass 
decreases from top to bottom.) The current-quark of 
perturbative QCD evolves into a constituent-quark 
as its momentum becomes smaller. The constituent- 
quark mass arises from a cloud of low-momentum 
gluons attaching themselves to the current-quark. 
This is DCSB: an essentially nonperturbative effect 
that generates a quark mass from nothing] namely, 
it occurs even in the chiral limit. 


provide a context, quantum electrodynamics, an Abelian gauge theory, has 12 672 diagrams 
at order a® in the computation of the electron’s anomalous magnetic moment 87|. Owing to 
its foundation in the non-Abelian group SU{3), the analogous perturbative computation of a 
quark’s anomalous chromomagnetic moment has many more diagrams at this order in the strong 
coupling. The number of diagrams represented by the self energy in Fig. [3] grows equally rapidly, 
i.e. combinatorially with the number of propagators and vertices used at a given order. Indeed, 
proceeding systematically, a computer will very quickly generate the first diagram in which the 
number of loops is so great that it is simply impossible to calculate in perturbation theory; 
impossible in the sense that we don’t yet have the mathematical capacity to solve the problem. 

Each of the diagrams which contributes to M(p^) in a weak-coupling expansion of Fig. [3] is 
multiplied by the current-quark mass, m. Plainly, any finite sum of diagrams must therefore 
vanish as m —)• 0. However, with infinitely many diagrams the situation might be very different: 
one has “0 x oo,” a product whose limiting value is contingent upon the cumulative magnitude 
of each term in the sum. Consider therefore the behaviour of M{p^) at large p^. QCD is 
asymptotically free [s^-l^. Hence, on this domain, each of the regularised loop diagrams must 
individually evaluate to a small number whose value depends on just how large is the coupling. It 
will not be surprising, therefore, to learn that for a monotonically-decreasing running-coupling, 
as{k‘^), there is a critical value of as( 0 ) above which the magnitude of the sum of infinitely many 
diagrams is sufficient to balance the linear decrease of m —)• 0 , so that the answer is nonzero and 
finite in this limit, viz., 


3 05 ( 0 ) I Vq; 5 ( 0 ) > 05 ( 0 ), Mq{p ) := lim M{p ; m) / 0. 

m —>-0 


(2.4) 


The internal consistency of QCD appears to guarantee that the limit is always finite. (As 
mentioned above, the case of Abelian theories is more complicated 2^ because they are not 
asymptotically free.) 


2.6. Challenge of Truncation 

The gap equation illustrates the features and flaws of each DSE. It is a nonlinear integral equation 
for the dressed-quark propagator and hence can yield much-needed nonperturbative information. 
However, the kernel involves the two-point function and the three-point function vl. The 
gap equation is therefore coupled to the DSEs satisfied by these functions, which in turn involve 
higher re-point functions. Hence the DSEs are a tower of coupled integral equations, with a 
tractable problem obtained only once a truncation scheme is specified. It is unsurprising that 










the best known truncation scheme is the weak coupling expansion, which reproduces every 
diagram in perturbation theory. This scheme is systematic and valuable in the analysis of 
large momentum transfer phenomena because QCD is asymptotically free but it precludes any 
possibility of obtaining nonperturbative information. 

Given the importance of DCSB in QCD, it is significant that the dressed-quark propagator 
features in the axial-vector WGTI, which expresses chiral symmetry and its breaking pattern: 

+ *["1/(0 +mg(C)]r{®(A:;P) = Sj^{k+)i'y5 + ij^Sg^ik-), (2.5) 

where /, g label quark flavours, P = pi +p 2 is the total-momentum entering the vertex and k is 
the relative-momentum between the amputated quark legs. To be explicit, k = {1 — rfjpi + r]P 2 , 
with g S [0,1], and hence k^ = pi = k + gP, k- = p 2 = k — {1 — g)P. In a Poincare covariant 
approach, such as presented by a proper use of DSEs, no observable can depend on g, i.e. 
the definition of the relative momentum. {N.B. Ref. [86| discusses the important differences 
encountered in treating flavourless pseudoscalar mesons.) 

In Eq. ([23]), T^ and Tg® are, respectively, the amputated axial-vector and pseudoscalar 
vertices. They are both obtained from an inhomogeneous Bethe-Salpeter equation (BSE), which 
is exemplified here using a textbook expression [91[: 

[r5^.{k; P)]tu = ^2[757/.]tn + r[S(g+)r 5 ^(<?; P)5((?_)],,R"(g, fc; P), (2.6) 

Jq 


in which K is the fully-amputated quark-antiquark scattering kernel, and the colour-, Dirac- and 
flavour-matrix structure of the elements in the equation is denoted by the indices r, s, t, u. N.B. 
By definition, K does not contain quark-antiquark to single gauge-boson annihilation diagrams, 
nor diagrams that become disconnected by cutting one quark and one antiquark line. In the 
parlance of Sec. II.41 K is two-particle irreducible (2PI). 

The WGTI, Eq. ()2.5p . entails that an intimate relation exists between the kernel in the 
gap equation and that in the BSE. (This is another example of the coupling between DSEs.) 
Therefore an understanding of chiral symmetry and its dynamical breaking can only be obtained 
with a truncation scheme that preserves this relation, and hence guarantees Eq. (12.51) without 
a fine-tuning of model-dependent parameters. Until 1995-1996 no one had a good idea about 
how to do this. Equations were truncated, sometimes with good phenomenological results and 
sometimes with poor results. Neither the successes nor the failures could be explained. 

That changed with Refs, il, 93|, which described a procedure that generates a Bethe-Salpeter 
equation from the kernel of any gap equation whose diagrammatic content is known. Its mere 
existence enabled the proof of exact nonperturbative results in QGD. (See, e.g.. Sec. V in 
Ref. 0 .) The scheme remains the most widely used today. Its leading-order term provides the 
rainbow-ladder (RL) truncation, which is accurate for ground-state vector- and isospin-nonzero- 
pseudoscalar-mesons 
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____, and the properties of ground-state octet and decouplet 

baryons [M-97] because corrections in these channels largely cancel, owing t o param eter-free 
preservation of the WGTIs. However, they do not cancel in other channels 9^101]. Hence 
studies based on the rainbow-ladder truncation, or low-ord er improv ements thereof, have usually 
provided poor results for scalar- and axial-vector-mesons 
states that are too low in comparison with other est imat es 
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sensitivity to model parameters for tensor-mesons [ll0| | and excited states |l0^ . llOOl . Illll . Ill2l | . 
In these circumstances one must conclude that physics important to these states is omit ted. 

Fortunately, a recently developed truncation scheme overcomes these difficulties ll 1,411 and is 
beginning to have a material impact. (An overview is presented in Sec. VI of Ref. |7ll].i This 
scheme is also symmetry preserving. Its additional strengths, however, are the capacities to 
work with an arbitrary dressed-quark gluon vertex and express DGSB nonperturbatively in the 
















































Bethe-Salpeter kernels. The new scheme has enabled a range of novel nonperturbative features of 
QCD to be demonstrated. For example, the existence of dressed-quark anomalous chromo- and 
electro-magnetic moments [HI and the key role they play in determining observable quantities 
j^; elucidation of t he c ausal connection between DCSB and the splitting between vector and 
axial-vector mesons 115] and the impact of this splitting on the baryon spectrum [^; and most 
recently, as explained in Sec. 12.81 below, a crucial step toward the ah initio prediction of hadron 
observables in continuum-QCD. 


2.7. Gluon Cannibalism 

I described DCSB and the dynamical generation of a running dressed-quark mass in Sec. 12.51 
Crucially, it is not just the propagation of quarks that is affecte d by strong interactions in QCD. 
The propagation of gluons, too, is described by a gap equation [llbl j: and its solution shows that 
gluons are cannibals: they are a particle species whose members become massive by eating each 
other! The associat ed g luon mass function, mg{k‘^), is monotonically decreasing with increasing 
and recent work [84] has established that 


mg{k‘^ =0) «0.5GeV. 


(2.7) 


The value of the mass-scale in Eq. (|2.7p is natural in the sense that it is commensurate with 
but larger than the value of the dressed light-quark mass function at far infrared momenta: 
M(0) ~ 0.3 GeV (see Fig.H]). Moreover, the mass term appears in the transverse part of the 
gluon propagator, hence gauge-invariance is not tampered with; and the mass function falls as 
l/k'^ for k'^ S> rng{0) (up to logarithmic corrections), so the gluon mass is invisible in perturbative 
applications of QCD: it has dropped to less-than 5% of it’s infrared value by k'^ = 4GeV^. 

Gauge boson cannibalism presents a new physics frontier within the Standard Model. 
Asymptotic freedom means that the ultraviolet behaviour of QCD is controllable. At the other 
extreme, dynamically generated masses for gluons and quarks entail that QCD creates its own 
infrared cutoffs. Together, these effects eliminate both the infrared and ultraviolet problems 
that typically plague quantum field theories and thereby make reasonable the hope that QCD 
is nonperturbatively well defined. 

The dynamical generation of gluon and quark masses provides a basi s for understanding 
the notion of a maximum wavelength for gluons and quarks in QCD 1171 ]. Indeed, given the 
magnitudes of the gluon and quark mass-scales, it is apparent that field modes with wavelengths 
A > ? ~ 2/mg{0) ~ 0.5fm decouple from the dynamics. They are screened in the sense described 
in Sec. [Til This is just one consequence of the appearance of a dynamically generated gluon 
mass-scale. 

There are many more. For example, the exceptionally light pion degree-of-freedom becomes 
dominant in QCD at those length-scales above which dressed-gluons and -quarks decouple 
from the theory owing to the large magnitudes of their dynamically generated masses. It is 
therefore conceivable that Gribov copies have no measurable impact on observables within 
the Standard Model because they affect only those gluonic modes whose wavelengths lie in 
the far infrared; and such modes are dynamically screened, by an exponential damping factor 
~ exp(—A/<j), so that their role in hadron physics is superseded by the dynamics of light-hadrons. 
This conjecture is consistent with the insensitivity to Gribov copies of the dressed-gluon and 
-quark tw o-point Schwinger functions observed in numerical simulations of QCD on fine lattices 
Ini,®. Another plausible conjecture is that dynamical generation of an infrared gluon mass- 
scale leads to saturation of the gluon parton distribution function at small Bjorken-x within 
hadrons. The possible emergence of this phenomenon stirs great scientific interest and curiosity. 
It is a key motivation in plans to construct an EIC that would be capable of producing a precise 
understanding of collective behaviour amongst gluons 
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Figure 5. Left panel - RGI running inte racti on strength computed via a combination of DSE- and 
lattice-QCD analyses, as explained in Ref. |llfil |. The function obtained with five different values of 
the renormalisation point is depicted in order to highlight that the result is RGI. The interaction is 
characterized by a value as(0) ~ 0.9 tt and the gluon mass-scale in Eq. (12.71) . Right panel ~ Gomparison 
between top-down results for the gauge-sector interaction (derived from the left-panel) with those obtained 
using the bottom-up approach based on hadron physics observables. Solid curve within grey band ~ 
top-down result for the RGI running interaction; and dashed curve within pale-green band ~ advanced 
bottom-up result obtained using the most sophisticated truncation of the matter sector DSEs - the 
DSE-DB kernel. The bands denote the domain of uncertainty in the determinations of the interaction. 


2.8. Continuum-QCD and ab initio predictions of hadron observables 

Within contemporary hadron physics there are two common methods for determining the 
momentum-dependence of the interaction between quarks: the top-down approach, which works 
toward an ab initio computation of the interaction via direct analysis of the gauge-sector gap 
equations; and the bottom-up scheme, which aims to infer the interaction by fitting data within 
a well-defined truncation of those equations in the matter sector that are relevant to bound- 
state properties. These two approaches have recently been united [s^ by a demonstration 
that the renormalisation-group-invariant (RGI) running-interaction predicted by contemporary 
analyses of QCD’s gauge sector coincides with that required in order to describe ground-state 
hadron observables using a nonperturbative truncation of QCD’s Dyson-Schwinger equations in 
the matter sector, i.e. the DCS B-improved (DSE-DB) kernel described briefly in Sec. 12.61 and 


elucidated in Refs. 71, 1131 - 115( |. 


The unification is illustrated in Fig. (5) the right panel presents a comparison between the 
top-down RGI interaction (solid-black curve within grey band, derived from the left panel) and 
the DB-truncation bottom-up interaction (green band containing dashed curve). Plainly, the 
interaction predicted by modern analyses of QCD’s gauge sector is in near precise agreement 
with that required for a veracious description of measurable hadron properties using the most 
sophisticated matter-sector gap and Bethe-Salpeter kernels available today. This is a remarkable 
result, given that there had previously been no serious attempt at communication between 
practitioners from the top-down and bottom-up hemispheres of continuum-QCD. It bridges a 
gap that had lain between nonperturbative continuum-QCD and the ab initio prediction of 
bound-state properties. 

It should be noted that if the realistic interaction depicted in Fig.[5] were employed as the 
seed for a RL-truncation study, it would fail completely because, inter alia, DCSB would be 
absent. We now know that a veracious description of DCSB and hence hadron properties in QCD 
requires a dressed-quark-gluon vertex. Constraining its form is a topic of great contemporary 
interest; and in this connection it cannot be emphasised too strongly that little of value today 
will be produced by any attempt at a term-by-term diagrammatic construction of this vertex. 






















2.9. Enigma of mass 

As noted in Sec. IM the pion is Nature’s lightest hadron. In fact, it is peculiarly light, with a 
mass just one-fifth of that which quantum mechanics would lead one to expect. This remarkable 
feature has its origin in DCSB. In quantum held theory the pion’s structure is described by 
a Bethe-Salpeter amplitude (here k is the relative momentum between the valence-quark and 
-antiquark constituents, and P is their total momentum): 

r^(A:;P) = 75 [iE^{k-,P) -b 7 • PF^{k-,P) + 7 • kG-^{k;P) - a^ykf,PyH.^{t, P)] , ( 2 . 8 ) 


which is simply related to an object that would be the pion’s Schrodinger wave function if a 
nonrelativistic limit were appropriate. I n QCD if, and only if, chiral symmetry is dynamically 
broken, then in the chiral limit (m = 0) 


UE^{k-,0) = B{k^ 


(2.9) 


where /^r is the pion’s leptonic decay constant, a directly measurable quantity that connects 
the strong and weak interactions, and the rhs is a scalar function in the dressed-quark prop¬ 
agator, Eq. (12.3p . This identity is miraculous. It means that the two-body problem is solved, 
almost completely, without lifting a hnger, once the solution to the one body problem is known. 
Eq. (|2.9p is a quark-level Goldberger-Treiman relation. It is also the most basic expression of 
Goldstone’s theorem in QGD, viz. 

Goldstone’s theorem is fundamentally an expression of equivalence between the one-body 
problem and the two-body problem in QCD’s colour-singlet pseudoscalar channel. 

Eq. (|2.9I) emphasises that Goldstone’s theorem has a pointwise expression in QGD; and, 
furthermore, that pion properties are an almost direct measure of the mass function depicted in 
Fig.[H Thus, enigmatically, properties of the (nearly-)massless pion are the cleanest expression 
of the mechanism that is responsible for almost all the visible mass in the Universe. Plainly, 
DCSB has a very deep and far-reaching impact on physics within the Standard Model. 


2.10. Nature’s Vacuum 

It will now be evident that DCSB is a crucial emergent feature of the Standard Model. It is 
very clearly expressed in the dressed-quark mass function of Fig.lH However, this perspective is 
relatively recent. DCSB was historically conflated with the existence of a spacetime-independent 
quark-antiquark condensate, (qq), that permeates the Universe. This notion was born with the 
introduction of QCD sum rule s as a theoretical artifice to estimate nonperturbative strong- 
interaction matrix elements P.fl^ and is typically tied to a belief, popularised in a Wikipedia 
entry, that the QCD vacuum is characterised by an infinite number of mass-dimensioned, 
spacetime-independent condensates, as illustrated in the left panel of Fig.[6l 

It is strange to imagine that in order to understand the properties of any single one of the 
hadrons supported by QCD, it is necessary to compute and explain the values of an infinite 
number of mass-dimensioned parameters, none of which is a property of any of the hadrons 
supported by QCD dynamics and all of which exist throughout the Universe, even in the absence 
of hadrons. This means, e.g. that the values of all these condensates in domains of the Universe 
which are not causally connected with ours nevertheless provide the keys to understanding 
the proton’s radius. This belief is too bizarre for me. Indeed, it is qualitatively equivalent 
to supposing that the Cooper pairs which explain superconductivity actually exist outside the 
superconducting metal instead of inside the lattice of po sitive ions that enables the electron- 


phonon interactions responsible for the phenomenon 122 l ]. 


There is an alternative, which can be stated succinctly as follows 121 - 125l |: 















Figure 6. Left panel - The QCD vacuum was historically imagined to be a “frothing sea,” with 
hadrons being merely bubbles of foam containing nothing but quarks and gluons interacting perturbatively 
throughout. Only near the bubble’s boundary did some so rt of stron g interaction occur, enforcing quark 
and gluon confinement. Right panel ~~ A modern view 12ll - ll25l| flips the picture completely. The 
space between hadrons is “empty,” except for the perturbative quantum fluctuations that characterise all 
interacting relativistic quantum field theories. The interior of hadrons, however, is extremely complex, 
with nonperturbative dynamics dominating in ~ 98% of the interior: the so-called condensates are 
spacetime-dependent and confined within the hadrons. 


If quark-hadron duality is a reality in QCD, then condensates, those quantities that were 
commonly viewed historically as constant empirical mass-scales that fill all spacetime, 
are instead wholly contained within hadrons, i.e. they are a property of hadrons 
themselves and expressed, e.g., in their Bethe-Salpeter or light-front wave functions. 

This standpoint is depicted in the right panel of Fig. [6] and canvassed fully in Sec. 4 of Ref. 

It presents the reasonable view that the understanding of hadrons requires that one explain 
what lies within those hadrons in contrast to the historical alternative, which suggested that 
hadrons could only be understood by explaining the properties of the vast spacetime domains 
that contain no hadrons at all. 

In order to provide a little mo re c ontext, it is perhaps worth noting that the Gell-Mann- 
Oakes-Renner (GMOR) relation 1261 ] is often cited as “proof” of the existence of a va cuum 


quark condensate. So let’s consider the published form of that identity, viz. Eq. (3.4) in Ref. |l26 |: 


ml = 




( 2 . 10 ) 


where ttItt is the pion’s mass and H^sh is that part of the hadronic Hamiltonian density which 
explicitly breaks chiral symmetry. It is crucial to observe that the operator expectation value in 
Eq. (|2.10l) is evaluated between pion states: it makes no reference to a vacuum quark condensate. 
Moreover, the virtual low-energy limit expressed in Eq. (|2.10l) is purely formal. It does not 
describe an achievable empirical situation. 

In terms of QCD quantities, Eq. (I2.10p entails 

~ 0, = m^^5^(0), Sl{0) =-{tt{P)\I{uu-I dd)\7r(P)) , (2.11) 

where with the running quark mass of the subscripted quark flavour, and 

5^(0) is the pion’s scalar form factor at zero momentum transfer, = 0. The rhs of Eq. (12.lip 




















is proportional to the pion a-term (see, e.g., Ref. |l27| |l. Consequently, using the connection 
between the cr-term and the Feynman-Hellmann theorem, Eq. (I2.10p is actually the statement 


~ 0 , ml = 


d 


dmi 




( 2 . 12 ) 


ud 


At this point one should recall the following mass formula for pseudo scala r mesons containing 
a valence /-quark and valence ( 7 -antiquark, which is exact in QCD [^, 12ll |: 


= {mj + m^g) , 
where, with the trace over colour and spinor indices, 

= (0|/g757^g/|Fo-) = ^2 trcD / i-i5l^lSf{k+)TH^_{k]P)Sg{k-), (2.14) 

J dk 

= -(0|/gi75g/|^^o-) = ^4 trcD / 'y 5 Sf{k+)rH^_{k]P)Sg{k-) (2.15) 

J dk 

are the meson’s pseudovector and pseudoscalar decay constants. 

Now, using Eq. (|2.13D . one obtains 


(2.13) 


5^(0) = 


d 




ml = 
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C P 


ud 
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(2.16) 


Equation (12.1611 is valid for any values of mu^, including the neighborhood of the chiral limit, 
wherein 


5^(0) = 


d 


dmid 


m 


C 

u 


ud 


p1_ 
/° ■ 


(2.17) 


The superscript “0” indicates that the quantity is computed in the chiral limit. With Eqs. (|2.11l) . 
()2.16p . (12.171) . one has shown that in the neighborhood of the chiral limit 


= -^id 


{qq)^° 




(2.18) 


This is a QCD derivation of the commonly recognised form of the GMOR relation. Neither the 
venerable tools of PCAC, i.e. exploiting the consequences of partial conservation of the axial 
current, nor soft-pion theorems were employed in analysing the rhs of Eqs. (12.101) . (12.lip . In 
addition, the derivation shows that in the chiral limit, the matrix element describing the pion- 
to-vacuum transition through the pseudoscalar vertex is equal to the expectation value of the 
explicit chiral symmetry breaking term in the QCD Lagrangian computed between pion states 
after normalisation by the pion’s decay cons tant . 

This recapitulation of the analysis in Ref. |l24l | emphasises anew that any connection between 
the pion mass and a “vacuum” quark condensate is purely a theoretical artifice. The true 
connection is that which one would expect, viz. the pion’s mass is a property of the pion, 
determined by the interactions between its constituents. One may further highlight the illogical 
nature of the vacuum connection by considering the expectation value 


{H{P)\qOq\H{P)), 


(2.19) 


where H{P) represents some hadron with total momentum P. If one chooses O = 7 ^, then 
all readers will accept that Eq. (I2.19p yields the electric charge of the hadron H. The choice 


































O = 757 ;i in Eq. (12.191) will yield the axial charge of the hadron H] O = will yield the tensor 
charge of the hadron H] and the choice 0 = 1 will yield the scalar charge of the hadron H. 
All these statements are true irrespective of the hadron involved, so that the scalar charge of 
the pion; i.e., 5|(0) = the “B-parameter” in chiral perturbation theory, is no more a 

property of the vacuum than is the pion’s electric charge. 

The change in understanding highlighted by shifting from the right to the left panel in 
Fig. [6] does not undermine the utility of the QCD sum rules approach to the estimation of 
hadron observables. Instead, it tames the condensates so that they return to being merely 
mass-dimensioned parameters in a useful computation scheme. Its implications are, however, 
significant and wide-ranging, especially if one imagines that the theory of gravity is understood 
well enough so t hat it may reliably be coupled to quantum field theory. Subscribers to this 
view argue that the energy-density of the Universe must receive a contribution from 

“vacuum condensates” and that the only possible covariant form for the energy of the (quantum) 
vacuum, viz. 

= ( 2 . 20 ) 


is mathem atica lly equivalent to the cosmological constant. From this perspective, the quantum 
vacuum is . a perfect fluid and precisely spatially uniform ...” so that “Vacuum energy 

is almost the perfect candidate for dark energy.” Now, if the ground state of QCD is really 
expressed in a nonzero spacetime-independent expectation value {qq)^ then the energy difference 
between the symmetric and broken phases is roughly Mqcd ~ 0.3 GeV, as indicated by Fig.Hl 
One obtains therefrom: 

^QCD ^ (2-21) 


i.e. the contribution from the QCD vacuum to the energy density associated with the 
cosmological constant exceeds the observed value by forty-six orders-of-magnitude. In fact, 
the discrepancy is far greater if the Hig gs vacuum expectation value is treated similarly. 


This mismatch has been called [13111 ] “... one of the gravest puzzles of theoretical physics.” 
However, it vanishes if one discards the notion that condensates have a physical existence, which 
is in dependent of the hadrons that express QCD’s asymptotically realisable degrees of freedom 
[T^ : namely, if one accepts that such condensates are merely mass-dimensioned parameters in 
one or another t heoretical computation and truncation scheme. This appears mandatory in a 
confining theory [T2^ - [T^ . a perspective one m ay em bed in a broader context by considering 


just what is observable in quantum field theory [I3ll ]: “.. .although individual quantum field 


theories have of course a good deal of content, quantum field theory itself has no content beyond 


analyticity, unitarity, cluster decomposition and symmetry.” If QCD is a confininj 
then the principle of cluster decomposition is only realised for colour singlet states 


_&eory, 
and 


all observable consequences of the theory, including its ground state, can be expressed via a 
hadronic basis. This is quark-hadron duality. Furthermore, if technicolour-like theories [1,S] are 
the correct scheme for explaining electroweak symmetry breaking, then the impact of the notion 
of in-hadron condensates is far greater still because it enables the Higgs vacuum expectation 
value to be understood and eliminated in the same manner. 


2.11. Parton Structure of Hadrons 

Since the advent of the parton model and the first deep inelastic scattering experiments there has 
been a d eterm ined effort to deduce the parton distribution functions (PDFs) of the most stable 
hadrons [l33( |. The behaviour of such distributions on the valence domain (Bjorken-x > 0.5) 
is of particular interest because this domain is definitive of hadrons, e.g. quark content on the 
valence domain is how one distinguishes between a neutron and a proton: a neutron possesses 
one valence u-quark plus two valence d-quarks whereas the proton possesses two valence u- 
quarks plus one valence d-quark. Indeed, all Poincare-invariant properties of a hadron: baryon 



















number, charge, total spin, etc., are d etermined by the PDFs which dominate on the valence 
domain. Moreover, via QCD evolution fT3i - [T^ . PDFs on the valence-quark domain determine 
backgrounds at the LHC. There are also other questions, e.g. rega rding flav our content of a 
hadron’s sea and whether that sea possesses an intrinsic component [1381 . Il39[ | . The answers to 
all these questions are essentially nonperturbative properties of QCD. 

Recognising the significance of the valence domain, a new generation of experiments, focused 
on Bjorken-x > 0.5, is planned at JLab, and under examination in connection with Drell-Yan 


studies at Fermilab and a possible EIC. Consideration is also being given to experiments aimed 
at measuring parton distribution functions in mesons at J-PARC. 

A concentration on such measurements requires theory to move beyond merely parametrising 
distribution functions and amplitudes. Computation within QCD-connected frameworks 
becomes critical because, without it, no amount of data will reveal anything about the theory 
underlying strong interaction phenomena. This is made clear by the example of the pion’s 
valence-quark PDF, u^{x), in connection with which a fa ilure of QCD was suggested foll owing 
a leading-order analysis of ttN Drell-Yan measurements [lia]. As explained in Ref. |l33| |. this 
confusion was fostered by the applica tion of a diverse range of models. On the other hand, a 
series of QCD-connected calculations |l4ll - 144| | subsequently established that the leading-order 
analysis was misleading, so that ul{x) may now be seen as a success for the unification of 
nonperturbative and perturbative studies in QCD. 

At this point it is worth noting that a meson’s Bethe-Salpeter wave function, which is obtained 
from its Bethe-Salpeter amplitude (see Eq. (|2.8p for the pion example) by reattaching the dressed- 
quark legs, viz. 

P) = Sfik+)rf<^ik- P)Sgik .), (2.22) 


is the quantum field theory analogue of the Schrodinger wave function that would describe the 
system if it were simply quantum mechanical, and whenever a nonrelativistic limit makes sense, 
the Bethe-Salpeter and Schrodinger wave functions become the same in that limit j9l|. For 
those desiring a probability interpretation of wave functions, this would be reassuring except, 
of course, for the fact that a nonrelativistic limit is never appropriate when solving continuum 
bound-state equations for composite systems containing the light u-, d- and s-quarks. 

To explain this comment, consider that the momentum-space wave function for a 
nonrelativistic quantum mechanical system, is a probability amplitude, such that 

|'0(p, t)P is a non-negative density which expresses the probability that the system is described 
by momenta p at a given equal-time instant t. Although the replacement of certainty in classical 
mechanics by probability in quantum mechanics was disturbing for some, the step to relativistic 
quantum field theory is still more confounding. As noted in connection with Eq. (|1.6I) . much of 
the additional difficulty owes to the loss of particle number conservation when this step is made. 
Two systems with equal energies need not have the same particle content because that is not 
conserved by Lorentz boosts and thus interpretation via probability densities is typically lost. 
To exemplify: a charge radius cannot generally be defined via the overlap of two wave functions 
because the initial and final states do not possess the same four-momentum and hence are not 
described by the same wave function. 

Such difficulties may be circumvented by formulating a theory on the light-front because the 
eigenfunctions of the light-front Hamiltonian are independent of the system’s four-momentum 
181-1201]. The light-front wave function of an interacting quantum system therefore provides a 


connection between dynamical properties of the underlying relativistic quantum field theory and 
notions familiar from nonrelativistic quantum mechanics. It can translate features arising purely 
through the infinitely-many-body nature of relativistic quantum field theory into images whose 
interpretation is seemingly more straightforward. Put simply; quantum mechanics-like wave 
functions can be defined in quantum field theory; quantum-mechanics-like expectation values 
can then also be defined and evaluated; and parton distributions may be obtained therefrom 



















Figure 7. Twist-two pion PDA computed 
using two very different DSE truncations at a 
scale C = 2 GeV. Dashed curve - rainbow-ladder 
(RL), the leading-order in a systematic, symmetry¬ 
preserving scheme 92, [^; and solid curve - the 


most sophisticated kernel that is currently available; 
namely, the DB kernel discussed in connection 
with Fig.[5l which incorporates nonperturbative 
effects generated by DCSB that are omitted in RL 
truncatio n and any stepwise improvement thereof 

11314115] , These results^ _ Me _ consistent with 

contemporary lattice-QCD [l46l - 149] . The dotted 
curve is = 6x( 1 — x), the result obtained 


in conformal QCD l4i 


via correlation functions evaluated at equal light-front time, namely within the initial surface 
x'^ = + x'^ = 0, and can thus be expressed directly in terms of ground-state light-front wave 

functions. Naturally, all that is only achieved if the light-front wave function can be calculated. 

The simplest object to compute is a hadron’s leading-twist two-particle parton distribution 
amplitude (PDA); and any framework that provides access to a hadron’s Poincare-covariant 
bound-state amplitude can also be employed to compute its PDAs. For example, the pion’s 
leading-twist two-particle PDA is given by the following projection of the pion’s Bethe-Salpeter 


wave function onto the light-front 145l ] 


/^V 9 ^(a:;C) = trcD-Z '2 / 5(n • g+ - x n • P) 757 • n S'(g+)r(g; P)5(g_), (2.23) 

J dq 


where: is a Poincare-invariant regularisation of the four-dimensional integral, with A 

the ultraviolet regularization mass-scale; Z 2 (CjA) is the quark wave-function renormalisation 
constant, with ( the renormalisation scale; n is a light-like four-vector, = 0; and P is the 
pion’s four-momentum, and n ■ P = —m.^-. 

Th e am plitude in Eq. (|2.23l) has been computed using two very different truncations of QCD’s 
DSEs (l45j |. with the result depicted in Fig. [3 Both kernels agree: compared with the asymptotic 
form, which is valid when ^ is extremely large, there is a marked broadening of ^nix‘, C)) which 
owes exclusively to DCSB. This causal connection may be claimed because the PDA is computed 
at a low renormalisation scale in the chiral limit, whereat the quark mass function owes entirely 
to DCSB via Eq. (|2.9p . Moreover, the dilation measures the rate at which a dressed-quark 
approaches the asymptotic bare-parton limit. It can be verihed empirically at JLabl2, e.g. 
in measurements of the pion’s electromagnetic form factor, the ratio of the proton’s electric 
and magnetic form factors, and the behaviour of the form factors that characterise transitions 
between the nucleon and its excited states. 

It is important to highlight that the computed PDAs in Pig. [3 are concave functions. Such 
pointwise behaviour contrasts markedly with the “humped” or “bimodal” distributions which 


have been favoured by some authors in p henomenological applications [ihfll ]. In this connection it 
must be understood, following Ref. that a double-humped form for the twist-two PDA lies 
within the class of distributions produced by a meson Bethe-Salpeter amplitude which may be 
characterised as vanishing at zero relative momentum, instead of peaking thereat. No ground- 
state pseu dosca l ar o r vect or meson Bethe-Salpeter equation solution exhibits corresponding 
behaviour [mi, imt. Thus, a bimodal distribution cannot be a pointwise-accurate 

representation of the PDA for a ground-state meson. On the other hand, if one is using such a 
























distribution in a practical phenomenological application for which only the lowest few moments 
are important, where the moments are defined via 


(x") = f dxx^(f{x), 

Jo 


(2.24) 


then some carefully-constrained bimodal distributions may serve to provide an approximation to 
the moments of a broad, concave PDA and can therefore provide useful information nevertheless. 

A question of more than thirty-years standing can be answered using Fig. [71 namely, when 
does (p^^{x) provide a good approximation to the pion PDA? Plainly, not at C = 2 GeV. The 
ERBL evolution eq uatio n l48l| describes the (^-evolution of (/9.n-(x;C); and applied to (p-n{x]C,) 


in Fig.lZl one finds [HfEi] that is a poor approximation to the true result even at 


C = 200 GeV. Thus at empirically accessible energy scales, the twist-two PDAs of ground-state 
hadrons are “squat and fat”. Evidence supporting this picture ha d Iq r ig b een accumulating 
I 52 I - 155 |: and the dilation is verified by simulations of lattice-QCD 1461 - 149l |. 


2.12. PDAs and Hard Exclusive Processes 

It is now possible to add flesh to the bones of Eq. (12.11) . In the theory of strong interactions, 
the cross-sections for many hard exclusive hadronic reactions can be expressed in terms of the 
PDAs of the hadrons involved. An example is the pion’s elastic electromagnetic form factor, for 
which the prediction can be stated succinctly 46l-l49l|: 


3Qo > Aqcd I Q^F-k{Q^) 


Q^>Qo 


167ras{Q‘^)f^'iV. 
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ipi 


OG = 


1 




1 


dx —(pnix ), 
lo x 


(2.25) 


where the running coupling is defined in Eq. (II.5p . As noted in Sec. 12.21 the value of Qo is not 
predicted by perturbative QGD. 

It was anticipated that JLab would verify this fundamental Standard Model prediction, 
Eq. (I2.25P : and in 2001, seven years after commencing operations, the facility provided the first 
high -precision pion electroproduction data for FV(Q ^) be tween values of 0.6 and 1.6 GeV^ 
156l |. In 2006 and 2007, t his d omain was revisited |l57l ] and a new result was reported, this 
time at = 2.45 GeV^ (l58| |. However, there was disappointment and surprise, with the 
collaboration stating that experiment is “still far from the transition to the region where the 
pion looks like a simple quark-antiq uark pair.” This data is represented in Eig.[8| by the filled 
circles and squares, drawn from Ref. [Hi, and the collaboration were comparing with Curve-D 
in the figure, which is the prediction obtained when ip^^{x), the asymptotic PDA appropriate 
to the conformal limit of QGD, is used in Eq. (12.251) . 

The confusion was compounded by the fact that JLab’s data confirmed th e be haviour 
predicted by a DSE-RL prediction for the pion form factor, computed in 2000 1 6,41 1. That 
prediction could not resolve the difficulty because it relied on brute numerical methods and 
hence could not produce a result at > 4GeV2 H In appearance, however, the shape of the 
prediction suggested to many that one might never see parton model scaling and QGD scaling 
violations at the momentum transfers ac cessib le to JLab, even after its upgrade. 

This conundrum was recently resolved [mil] . Using a refinement of known methods 
also employed in the successful analysis of (/? 7 r(x;C), described in Sec. 12.Ill a reliable prediction 
of is now available on the entire domain of spacelike Q^. This is Gurve-A in Fig.jSJ 

Moreover, the analysis enables correlation of that result with Eq. (|2.25l) . using the modern 
PDA computed in precisely the same framework, which is Gurve-G in the figure. This leading- 
order, leading-twist QGD prediction, obtained with a pion valence-quark PDA evaluated at a 


® Notably, modern IQCD computations are restricted to the domain < 3GeV^ [l64| . 













































Figure 8. Solid curve (A) - 

Theoretical prediction |16Cll| : dotted curve (B) 

- monopole form fitted to data; dot-dot-dashed 
curve (C) - perturbative QCD (pQCD) prediction 
computed with the modern, dilated pion PDA 
described in Sec. l2.11l and dot-dot-dashed curve (D) 

- pQCD prediction computed with the conformal- 
limit PDA, which had previously been used to 
guide expectations for the asymptotic behaviour of 

The fil led-c ircles and -squares represent 
existing JLab data [l59l| ; and the filled diamonds and 
triangle, whose normalisation is arbitrary, indicate 
the projected Qf-reach and accuracy of forthcoming 


experiments 161 


scale appropriate to the experiment, underestimates the full DSE-RL computation by merely 
an approximately uniform 15% on the domain depicted. The small mismatch is explained 
by a combination of higher-order, higher-twist corrections to Eq. (12.2511 in pQCD on the one 
hand and, on the other hand, shortcomings in the rainbow-ladder truncation (see Sec. l2.8p . 
which predicts the correct power-law behaviour for the form factor but not precisely the right 
anomalous dimension (exponent on the logarithm) in the strong-coupling calculation. Hence, 
disappointment is now transformed into optimism because the comparison of Curves A and C 
in Eig. [5] predicts that the upgraded JLab facility will reveal a maximum at ~ 6GeV^ and 
an experiment at = 9 GeV^ will see a clear sign of parton model scaling for the first time in 
a hadron elastic form factor. 

The implications of these results are far-reaching. They indicate that hadro-particle physics 
is on the cusp of verifying the theory of factorisation in hard exclusive processes, with dominance 
of hard contributions to on > 8GeV^. Notwithstanding that, the normalisation of 

this form factor is fixed by a pion wave-function whose dilation with respect to = 6x{l — x) 

is a definitive signature of DCSB; and this affords the opportunity for an empirical measurement 
of the strength of DCSB, i.e. the power behind the origin of visible mass. This series of outcomes 
will be important pages in a book on the Standard Model, in which the first lines were written 
forty years ago; and they pave the way for a dramatic reassessment of pictures of proton and 
neutron structure, which is already well underway, as we shall see in Sec.O 


2.13. Light-front Wave Function 

A modern perspective on the pion’s PDA was presented in Sec. 12 .Ill The more general quantity 
is the pion’s light-front wave function: 


f^(fT,{x,k‘i) = Z2tTcD j ■ q+- xn-■nS{q+)r{q;P)S{q-), (2.26) 

from which one may readily recover the PDA by integrating: f dqidq 2 , with a careful redefinition 
of the renormalisation constant and regularisation scale, as is plain by comparison wit h 
Eq. (I2.2.sp . If one proceeds from Eq. (12.2611 using precisely the analytical framework in Ref. 145l | , 
then one will obtain a factorised result, viz. 


ip^{x, kl) = ipl{x) X ipl{k \), 


(2.27) 


an outcome which cannot be true in any nontrivial theory. 






















To expose how this factorisation arises, consider the following illustrative, algebraic 
representations for the propagators and amplitudes: 


S{p) 

= [-Z7 •p + M]AAf(p^), 

(2.28a) 

Pu{z) 

1 r(z^ + 3/2) ^2\u 

A r(i/ + i) ^ ^ ’ 

(2.28b) 




r^(A:;g) 

- n5 f / dzp{z)AM{k+^), 

J-K J-1 

(2.28c) 


where Am{s) = l/[s + k±z = /c — (1 ^ z)qj 2 , and I have set = 0 and used a momentum 
partitioning parameter ij = 0. Inserting these expressions into Eq. (I2.26p . one obtains an 
integrand denominator that is a product of fe-quadratic forms, each raised to some power. The 
denominator product can be combined into a single quadratic form, raised to a unique power, 
by using a Feynman parametrisation, so that Eq. (I2.26p yields 


(x^ikD) 


dz p,y{z)u{i/ + 
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X Z2trcD 
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^ d'^kii 




n ■ k 


dxidx2 xi{l — xi) 


v-l 




(27r)2 \n ■ P J [{kd — f{xi,X 2 , z)Py + k'j_ + ^ 


(2.29) 


where k^^ = {ks,ki) and f{xi,X 2 ,z) is some simple, algebraic function of its arguments. At this 
point one would shift variables to obtain 

iq + q + MT*'- 

X / dz puiz)i'{i'+ 1) 


f 


dxidx 2 Xi{l - xi)'-' ^[/(xi,X2,2)]”* , (2.30) 


and the factorisation is now apparent: the rhs of the first line in this equation is independent of 
the second line, which contains all nontrivial information about the moments. 

In order to learn how one might correct this, return to the expression for the PDA’s moments: 





^Tr{x, A:_i_), 


(2.31) 


namely. 


(x^) ~ J dz pv{z)v{v + 1) J dxidx 2 Xi{l — xi^ ^ 
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= ■^2tl'CD 
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^ d^k 


(27r)4 [p + M2]2+i^ 

J dz Py{z)v{v+ 1) J dxidx2Xi{l - xiy~^[f{xi,X2,z)Y 


(2.32) 


(2.33) 



















Now compare this with a slightly modified expression: 


J dz p„{z){l - + 1 ) j dxidx 2 Xi{l - xiY ^ 

d'^k f n ■ k 
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X Z2trcD 


M 


= ■^2trCD 
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(27rY \n ■ P J [{k — f{xi,X 2 ,z)PY + {1 — 

^ d^k 

[fc2 + (1 - 22)M2]2+^ 

X J dz p^{z){l - z'^Y u{p + 1) J dxidx 2 Xi{l - xiY~^[f{xi,X 2 ,z)]^. 

A simple change of variables: —)• = (1 — converts this expression into 

{x'^)m ~ J dz pYz)Y - z'^Y + Y J dxidx2Xi{l - xiY~^[f{xi,X2,z)Y 


(2.34) 


(2.35) 


xZ2trcD 


= xZ2trcD 


d^k 


2\2 


{1-zY 




(27r)4'" [k‘^{l-zY + {l-zYM^Y+^ 
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(2.36) 


(27r)4 [k^ + A/2]2+^ 

J dz Pu{z) v{v + 1) J dxidx 2 Xi{l - xiY~^[f{xi,X 2 ,z)]"^ (2.37) 

= (x™). (2.38) 

Plainly, so far as a computation of the PDA’s moments is concerned, Eq. (12.341) is actually the 
same as Eq. (j2.32p . However, a difference appears if one reverts to the light-front wave function 
in Eq. ()2.26p and begins with Eq. (I2.34p . viz. 
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In this way one arrives at 
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2u 


[M2 + A:2/(1-z2)]1+^ 


(2.43) 

(2.44) 

(2.45) 


and now, very clearly, {x^{k‘j_)) has not a factorised: the dz integral is modulated by the 
k\/{l — z'^) part of the denominator. Consequently, (pT^{x,k'\) will not be expressible as a 
product of the type in Eq. (|2.27l) . 

In order to provide a practical illustration, let’s work with a simple expression for the pion’s 
Bethe-Salpeter wave function: 

xAk\ P) = S{k)T{k- P)S{k -P)= i 75 S^{k, k-P)- ^757 • PP^{k, k - P). (2.46) 

JTT 

I have omitted the Gnik', P) and P^ik] P) terms for the sake of simplicity. That won’t alter the 
gist of the following analysis. In terms of Xm the moments of the pion’s PDA are: 

fn (n ■ P)™+i (x”") = trcD-Z '2 f (n • /c)™ 757 • n Xnik; P) (2.47) 

Jdk 

= ^AN,Z 2 [ {n-kr{n-P)P^{k,k-P), (2.48) 

In J dk 


which is a very simple expression, that further reduces to 

f-A 


fl {n • pr {xn = 4 iVeZ 2 [ (n • krP^{k, k-P) 

Jdk 


(2.49) 


Focusing on the PDA’s zeroth moment, one has 


i-A 

/2 = 4IVe^2 / PAk,k-P). 

Jdk 


(2.50) 


One knows that in Landau-gauge Z 2 (C,A) behaves as [1/ In In A] for Aqcd/A ~ 0 |l2ll |. This 
means that a nonzero value for is possible if, and only if, 


Pnik, k — P) 


fe^^AqcD 
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A:2 lnA:2/Aq 


(2.51) 


QCD 


It follows that if one wishes to make a model of Pn{k, k — P) for use in algebraic analysis, then 
a QCD-like theory will have 


PTr{k,k-P) = J dzpi{z) 


AP 


[{k - ?P)2 + M2]2 ’ 


(2.52) 













where = (1 — 2 )/ 2 . Inserting this into Eq. (j2.5np . one hnds 


= ANcZ2 
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,,[fe 2 +M 2]2 7_ 


J dzpi{z) 


= 4iV,Z2 


m 2 


',,[fc 2 + M 2 ] 2 - 


(2.53) 


With Z 2 —>• 1, this is readily re cogn ised as the expression for /2 that is obtained in models of 
the Nambu-Jona-Lasinio type 167 | when an 0(4)-symmetric regularisation scheme is nsedP^ 
If one restores the logarithms, then it’s a QCD-like expression. 

Consider now the expression for all PDA moments: 


k —/* 7 / \ 

= J dzpi{z)> 
It follows immediately that 


M 2 
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6 

{m + 2)(m + 3) 


<f{x) = = 6 x(l — x ), 


(2.54) 

(2.55) 

(2.56) 


i. e. Eq. (j2.46l) is sufficient to reproduce the result hrst explicated in Ref. |l45l | . 

It is now time to return to the light-front wave function, defined via its moments, as before; 
but corrected as described above. Using 
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M 2 
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(2.58) 

(2.59) 

(2.60) 


One may now analyse {x'^{kj_)) via a twist expansion, viz. an expansion in 1 /fe 


2 . 


1 N P 

= -J 2 ^ j ^dzp,{z),^‘ 


m 2 , „ M^ 

■p (1 “ 2 ; )-^ + 


fe^ 


(2.61) 


Eocusing hrst on the leading-twist coefficient - the factor accompanying l/k\, one has 

(x"*(fei))( 2 ) = 


As this is being written, I have learnt that Prof. Yoichiro Nambu died on 5 July 2015, at the age of 94. Winner 
of the Nobel Prize in 2008 for discovering the mechanism of DCSB, he was a great physicist; but also a fine 
person. I have learnt from people who knew him well that Prof. Nambu was always deeply interested in exotic 
phenomena and intent on gathering related information from a broad range of sources. Such habits distinguish 
the best in this discipline. 
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Figure 9. Left panel - Profile of light-front wave function in Eq. ()2.7np . k^), computed 

with M = 0.4 GeV. This expression provides an “asymptotic” pion wave function, in a sense 
similar to that in which Eqs. (I2.28p p rodu ce an asymptotic PDA. Right panel - Profile of the 
model wave function described in Ref. the additional curvature is typical of wave functions 

relevant for the description of phenomena explored at a typical hadronic scale, e.g. 1 GeV. 


from which it is plain that 
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/' 27r2 

The subleading twist coefficient ~ the factor accompanying 1 /k\_ - is 
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and hence 




fn 3vr2 

One may readily establish that this pattern continues. 

At the other extreme, in the neighbourhood k'^ ~ 0, one finds: 

(x^(k‘^ = 0)) =f dzo](z)c^ ^ _ 1 Ve 3 1 

^ f^27r^ 1-^2 /2 27r2 2 1 + m’ 


(2.64) 

(2.65) 

( 2 . 66 ) 

(2.67) 

( 2 . 68 ) 


namely, 

ip^{x,k‘j_ =0) = - x), (2.69) 

where 77(1 — x) is the Heaviside step-function. This is a natural result, at least qualitatively: 
it means that on /c2 ~ Q the valence-quarks are free to roam with (almost) no bound: each 
X G [0,1] has equal probability. N.B. In QGD-like theories, this translates into a very broad 
distribution, that is nonetheless zero at x = 0,1. 

With this information in hand, one is able to reconstruct the expression for the pion light-front 
wave function produced by Eqs. (I2.39p . (I2.46p . (I2.57p : 


f^ip^{x,kl) 


_9_]_ 

47r2 k‘^. 

1 +-^ 


(2.70) 


4M2x( 1 — x) 
























which is depicted in Fig.[9l This expression exposes the natural form for the correlated 
dependence of the two-body wave function on {x,k\), i.e. it depends on the two-body light- 
front kinetic energy: k\/[x{l — x)]. 


3. Baryons as they really are 
3.1. Grand Unification 

There are many novel challenges in the Standard Model to which one might now turn. For 
example, a QCD-connected prediction of the spectrum of hybrid and exotic mesons would be 
very valuable given that a flagship effort at JLab 12 is the search for these states. Equally, 
or perhaps more pressing is the need to address the three valence-quark bound-state problem 
in QCD with the same level of sophistication that is now available for mesons, with the goal 
being to correlate the properties of meson and baryon ground- and excited-states within a single, 
symmetry-preserving framework. Here, symmetry-preserving means that the analysis respects 
Poincare covariance and satishes the relevant WGTIs. Constituent-quark models have hitherto 
been the most widely applied spectroscopic tools; and whilst their weaknesses are emphasised by 
critics and acknowledged by proponents, they are of continuing value because there is nothing 
better that is yet providing a bigger picture. Nevertheless, they possess no connection with 
quantum held theory and therefore no connection with QCD; and they are not symmetry¬ 
preserving and therefore cannot voraciously connect meson and baryon properties. 

A truly comprehensive approach to QCD will provide a unihed explanation of both mesons 
and baryons. I have explained that DCSB is a keystone of the Standard Model, which is 
evident in the momentum-dependence of the dressed-quark mass function - Fig.HJ it is just as 
important to baryons as it is to mesons. The DSEs furnish the only extant framework that can 
simultaneously and transparently connect both meson and baryon observables with this basic 
feature of QCD, having provided, e.g. a direct correlation of meson and baryon properties via a 
single interaction kernel, which preserves QCD’s one-loop renor malisation group behaviour and 
can systematically be improved. This is evident in Refs. 


3.2. Borromean Analogy 

Let us focus initially on the proton, which is the core of the hydrogen atom, lies at the heart 
of every nucleus, and has never been observed to decay. It is nevertheless a composite object, 
whose properties and interactions are determined by its valence-quar k co ntent: u u + d, i.e. 
two up (u) quarks and one down (d) quark. So far as is now known [l75l |. bound-states seeded 
by two valence-quarks do not exist; and the only two-body composites are those associated with 
a valence-quark and -antiquark, i.e. mesons. These features are supposed to derive from colour 
confinement, whose complexities are discussed in Sec. 12.41 

Such observations lead one to a position from which the proton may be viewed as a Bor¬ 
romean bound-state j^, viz. a system constituted from three bodies, no two of which can 
combine to produce an independent, asymptotic two-body bound-state. In QCD the complete 
picture of the proton is more complicated, owing, in large part, to the loss of particle number 
conservation in quantum held theory and the concomitant frame- and scale-dependence of any 
Fock space expansion of the proton’s wave function which is described briehy in Sec. II.21 

Notwithstanding that, the Borromean analogy provides an instructive perspective from which 
to consider both quantum mechanical models and continuum treatments of the nucleon bound- 
state problem in QCD. It poses a crucial question: 


Whence binding between the valence quarks in the proton, i.e. what holds the proton 
together? 


In numerical simulations of IQCD that use static sources to represent the proton ’s va l ence - 
quarks, a “Y-junction” hux-tube picture of nucleon structure is produced, e.g. Ref. [T^.[r^. 













This might be viewed as originating in the three-gluon vertex, which signals the non-Abelian 
character of QCD and is the source of asymptotic freedom [881490 1 as discussed in connection 
with Fig.[T] Such results and notions would suggest a key role for the three-gluon vertex in 
nucleon structure if they were equally valid in real-world QCD wherein light dynamical quarks 
are ubiquitous. However, as we saw in Sec. 12.41 they are not; and so a different explanation of 
binding within the nucleon must be found. 

DCSB has numerous corollaries that are crucial in determining the observable features of the 
Standard Model, some of which are detailed above. Another particularly important consequence 
is less well known. Namely, any interaction capable of creating pseudo-Goldstone modes as 
bound-states of a light dressed-quark and -antiquark, and reproducing the measured value of 
their leptonic decay constants, will necessarily also generate strong colour-antitriplet correlations 
between any two dressed quarks contained within a nucleon. Although a rigorous proof within 
QCD cannot be claimed, this assertion is based upon an accumulated body of evidence, gathered 
in two decades of st udying two- and three- body bound-state problems in hadron physics, e.g. 
Refs. [13, [Hi, [13, m, [Iffl, Eoll, [TtO, [iTi-EIi] . No realistic counter examples are known; an d the 
existence of such diquark correlations is also supported by simulations of IQCD isd . 187l ]. 

The properties of diquark correlations have been charted. Most importantly, diquarks are 
confined. How ever, this is not true if the RL truncation is used to define the associated scattering 
problem 182l |. Corrections to that simplest symmetry-preserving approximation are critical in 
quark-quark channels: they eliminate bo und-state poles from the quark-quark scattering matrix 
but preserve the strong correlations [13, HH. 

Additionally, owing to properties of charge-conjuga tion , a diquark with spin-parity 


may 


be viewed as a partner to the analogous J ^ meson [Tt^ . R follows that scalar, isospin-zero 
and pseudovector, isospin-one diquark correlations are the strongest; and whilst no pole-mass 
exists, the following mass-scales, which express the strength and range of the correlation and 
are each bounded below by the partnered meson’s mass, may be associated with these diquarks 

[l^ . f^ . 


m 


Mo+ ~ ~ ’ rn{uu}^+ ~ 0.9 - 1.1 GeV , 


(3.1) 


with in the isospin symmetric limit. Realistic diquark correlations 

are also soft. They possess an e lectromagnetic size that is bounded below by that of the analogous 
mesonic system, viz. 188 . 189( |: 


r[ud]o+ ^ ^ rp , (3.2) 

with . As with mesons, these scales are set by that associated with DCSB. 

It is useful to remark that diquarks are colour-singlets in a dynamical theory based on 
SU(2)-colour. They would thus exist as asymptotic states and form mass-degenerate multiplets 
with mesons composed from like -flavoured quarks. (These properties are a manifestation of 
Pauli-Giirsey symmetry .) Gonsequently, the [ud]Q+ diquark would be massless in the 

presence of DGSB, matching the pion, and the {ud}i+ diquark would be degenerate with the 
theory’s p-meson. Such identities are lost in changing the gauge group to SU(3)-colour; but 
clear and instructive similarities between mesons and diquarks nevertheless remain. 

3.3. Diquarks in the nucleon 

It should now be apparent that the bulk of QGD’s particular features and nonperturbative 
phenomena can be traced to the evolution of the strong running coupling. Its unique 
characteristics are primarily determined by the three-gluon vertex; the four-gluon vertex does 
not contribute dynamically at leading order in perturbative analyses of matrix elements; and 
nonperturbative continuum analyses of QGD’s gauge sector indicate that satisfactory agreement 

















































Figure 10. Poincare covariant Faddeev equation. 

is the Faddeev amplitude for a baryon of total 
momentum P = pq+pd, where Pq^d are, respectively, 
the momenta of the quark and diquark within 
the bound-state. The shaded area demarcates the 
Faddeev equation kernel: single line, dressed-quark 
propagator; F, diquark correlation amplitude; and 
double line, diquark propagator. 


with gluon propagator results from IQCD simulations is typicall y o b tain e d withou t reference to 
dynamical contributions from the four-gluon yertex, e.g. Refs.]?^. 116 . 192 - 199l |. The three- 
gluon vertex is therefore the dominant factor in producing the class of RGI running interactions 
that have provided both successful descriptions and predictions of many hadron observables 
It is this class of interactions that generates the strong attraction between 


71 


[h 

two quarks which produces tight diquark correlations in analyses of the three valence-quark 
scattering problem. 

The existence of tight diquark correlations considerably simplifies analyses of the three 
valence-quark scattering problem and hence baryo n bo und states because it reduces that task 
to solving a Poincare covariant Faddeev equation flT^ . depicted in Fig.lTOl The three gluon 
vertex is not explicitly part of the bound-state kernel in this picture of the nucleon. Instead, 
one capitalises on the fact that phase-space factors materially enhance two-body interactions 
over n > 3-body interactions and exploits the dominant role played by diquark correlations in 
the two-body subsystems. Then, whilst an explicit three-body term might affect hne details of 
baryon structure, the dominant effect of non-Abelian multi-gluon vertices is expressed in the 
formation of diquark correlations. Such a nucleon is then a compound system whose properties 
and interactions are primarily determined by the quark-bdiquark structure evident in Fig. 1101 

It is important to highlight that both scalar-isoscalar and pseudovector-isotriplet diquark 
correlations feature within a nucleon. Any study that neglects pseudovector diquarks is 
unrealistic because no self-consistent solution of the Faddeev equation in Fig. (TU] can produce a 
nucleon constructed solely from a sc alar diquark, e.g. pseudovector diquarks typically provide 
roughly 150 MeV of attraction 


183|. 


The relative probability of scalar versus pseudovector 
diquarks in a nucleon is a dynamical stat eme nt. Realistic computations predict a scalar diquark 
roximately 60% [ 


64l. Il85l |. This prediction can be tested by contemporary 


strength of apj 
experiments |65| |. 

The quark-kdiquark structure of the nucleon is elucidated in Fig. 1111 which provides a 
representation of the leading component of the nucleon’s Faddeev amplitude: with the notation 
of Ref. [6^, si(|p|,cos0), computed using the Faddeev kernel described therein. This function 
describes a piece of the quark-l-scalar-diquark relative momentum correlation. Notably, in this 
solution of a realistic Faddeev equation there is strong variation with respect to both arguments. 
Support is concentrated in the forward direction, cos 9 > 0, so that alignment of p and P 
is favoured; and the amplitude peaks at {\p\ ~ M]\f/6,cos6 = 1), whereat pq « P/2 Ri p^ 
and hence the natural relative momentum is zero. In the antiparallel direction, cos 0 < 0, 
support is concentrated at \p\ = 0, i.e. pq « F/S, pd ~ 2P/3. A realistic nucleon amplitude is 
evidently a complicated function; and significant structure is lost if simple interactions and/or 
truncations are employed in building the Faddeev kernel, e.g. extant treatments of a momentum- 
independent quark-quark interaction - a contact interaction - produce a Faddeev amplitude that 
is also momentum independent (^ . 631, a result exposed by Fig. [11] as unrealistic for any probe 
sensitive to the nucleon interior. 

A nucleon (and kindred baryons) described by Fig.lTOjis a Borromean bound-state, the binding 
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Figure 11. Representation of the dominant 
piece in the nucleon’s eight-component Poincare- 
covariant Faddeev amplitude: si(|p|,cos0). In 
the nucleon rest frame, this term describes that 
piece of the quark-diquark relative momentum 
correlation which possesses zero intrinsic quark- 
diquark orbital angular momentum, i.e. L = 0 
before the propagator lines are reattached to form 
the Faddeev wave function. Referring to Fig. (TUI 
p = P/3 — pq and cos9 = p ■ Pj^JjFP^. (The 
amplitude is normalised such that its Pq Chebyshev 
moment is unity at \p\ =0.) 


within which has two contributions. One part is expressed in the formation of tight diquark 
correlations; but that is augmented by attraction generated by the quark exchange depicted in 
the shaded area of Fig.lTOl This exchange ensures that diquark correlations within the nucleon 
are fully dynamical: no quark holds a special place because each one participates in all diquarks 
to the fullest extent allowed by its quantum numbers. The continual rearrangement of the 
quarks guarantees, inter alia, that the nucleon’s dressed-quark wave function complies with 
Pauli statistics. 

It is impossible to overstate the importance of appreciating that these fully dynamical diquark 
correlatio ns are vas tly different from the static, pointlike “diquarks” which featured in early 
attempts 20 Cll . 20 ll | to unde rstand the baryon spectrum and to explain the so-called missing 
resonance problem 202H204i | . Modern diquarks are soft, Eq. (|3.2I) : and enforce certain distinct 
interaction patterns for the singly- and doubly-represented valence-quarks within the proton. 
On the other hand, the number of states in the spectrum of baryons obtained from the Faddeev 
equation in Fig.[TU] (205j | is similar to that found i n the three-constituent quark model, just as it 
is in today’s IQCD calculations of this spectrum (206l |. 


3.4- Contact interaction 

It is only very recently that numerical algorithms have been developed that both enable one 
to use sophisticated DSE kern e ls in order to gain direct access to the large-Q^ behaviour of 
hadron form factors (l45l . 16(i . 207 ] and promise to pave the way for their use in treating 
the phenomena of deep inelastic scattering. Absent those algorithms, a confining, symmetry¬ 
preserving DSE treatme nt o f a ve c tor 0 vec tor contact interaction has proved useful in a variety 
of contexts (gH, Il8,3l . Il89 . 1208 - 21,31 1. This collection of work is dwarfe d, however, by 


the widespread application of models of the Nambu-Jona-Lasinio (NJL) type (l67l . 214l - 217 |. 
It must be remarked that whilst there is a ground-level similarity, these two approaches are 
different. NJL-model practitioners allow themselves the freedom of tuning the relative strength 
of individual four-fermion interaction terms in a model Lagrangian in order to fit a body of data 
with the aim of providing a phenomenology of hadron physics that can serve both to elucidate 
correlations between observables and as a beacon that sheds light on novel explanations for 
unexpected phenomena. The DSE approach, on the other hand, treats the contact interaction 
as a representation of the gluon’s two-point Schwinger function. At the outset, therefore, the 
DSE approach fixes the number, t ype and relative strength of the four-fermion interaction terms, 
illustrated elsewhere Hii,!!!, and then proceeds to identify and highlight those observables 


as 


that can distinguish between different choices for the momentum-dependence of the DSE kernels. 
There is true predictive and discriminative power in such analyses and hence this is the class of 































































189l . 208, 200( 1 showed that a confining, symmetry-preserving 


applications that I will illustrate. 

The earliest studies in this class 
treatment of a vector 0 vector contact interaction provides results for vr- and p-meson observables 
that are not realistically distinguishable from those obtained with the best renormalisation- 
group-improved one-gluon exchange interactions so long as the momentum of the probe involved, 
Q, does not exceed the zero-momentum value of the dressed-quark mass, i.e. on 
Furthermo re, w hilst a contact interaction typically produces form factors that are too hard 
[hll, [H, [ 93 , , interpreted judiciously, even these results can be used to draw valuable insights, 

e.g. concerning the relationships between different hadrons. Importantly, too, studies employing 
a symmetry-preserving regularisation of the contact interaction serve as a useful surrogate, 
enabling the exploration of domains which analyses using interactions that more closely resemble 
those of QCD are as yet unable to enter. They are therefore useful at present in attempts to 
use data as a tool for charting the nature of the quark-quark interaction at long-range, i.e. for 
identifying clear signals in observables for the running of couplings and masses in QCD that was 
elucidated in Secs. 12.5112.7112.81 

Some of the more interes ting and robust studies are contact-interaction analyses which unify 
meson and baryon spectra 1831 . 205| |. Given that I have already described the gap and Bethe- 


Salpeter equations, which are crucial in computing the meson spectrum, at this point it is worth 
illustrating a Faddeev equati on; a nd the simplest example is that obtained for the A-baryon 


using the contact-interaction 1831 ] : 


_ ogj ^991+ f ^ + (1 " a?rn\){am/^ + M) 

M J Jq “ + cTA{a,M,mgg^^,mAW 

where C™{z) = C^i{z)/z, C™{z) = —z{d/dz)C™{z), — F(— 1 , with 

T(a,y) being the incomplete gamma-function, and 

(TA(a, M, mqq^j _, ttia) = (1 — a) -|- a rriqq^j^ — a{l — a) m\. (3.5) 

Equation (13.41) is an eigenvalue problem whose solution yields the mass for the dressed-quark- 
core of the A-resonance. The primary parametric inputs to this equation are the: dressed-quark 
mass, M; axial-vector diquark mass, mqq^q_, and correlation strength, Eqq^q_; and the contact- 
interaction regularisation parameters, Tir, Tuv The former are computable once the latter are 
specified. A final parameter is qa, which appears as the result of a drastic simplification, viz. in 
the Faddeev equation for a baryon of type B = N, A, the quark exchanged between the diquarks 
in the shaded region of Fig.[T 0 ]is represented as 


(3.3) 

(3.4) 


(3.6) 

This is a variant of the so-called “static approximation,” which itself was in trod uced in Ref. |22[ll ] 
and has subsequently been used in studies of a range of nucleon properties mu. In combination 
with diquark correlations generated by the contact-interaction, whose Bethe-Salpeter amplitudes 
are momentum-independent, Eq. (13.6p generates Eaddeev equation kernels which themselves are 
momentum-independent. The dramatic simplifications which this produces are the merit of 
Eq. (|3.6I) . The analogous equation for the nucleon is a five-dimensional algebraic eigenvalue 
problem because the nucleon also contains scalar diquark correlations. 

In order to proceed it is important to highlight just what is missing in all DSE kernels that 
are currently available. Namely, even the best kernels are built only from dressed-gluons and 
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Figure 12. Left panel: Masses of m embers of the baryon octet: Circles - computed masses [2 Of 
and diamonds - empirical masses 175l |. The horizontal axis lists a particle name with a subscript, and 


reading left-to-right the labels indicate ground-state, radial-excitation, ground-state’s parity partner and 
radial excitation of that state. Right panel: Analogous plot for the decuplet masses. 


-quarks. They omit long-range interactions. However, QCD produces a very potent long-range 
interaction, viz. that associated with the pion and other light mesons, without which no nuclei 
would be bound and we wouldnt be here. Fortunately, since the real character of contemporary 
DSE kernels is understood, one knows that such resonant “meson-cloud” effects can be added 
without over-counting. Uncovering just how that can be accomplished in a symmetry -preservin g 
manner is where questions still lie. Only rudimentary attempts exist {e.g. Refs. (2221 - 22^ ). 
Notwithstanding that, meson-cloud contributions must be thoughtfully considered before any 
comparison may be made between contemporary computations and the real world. In this 
connection it is now judged that whilst resonant contributions to truncated DSE kernels may 
have a material impact on the masses of the nucleon and A-baryon separately, the modification 
of each is approximately the same, so that the mass difference, 5m = — m^, is la rgely 

unaffected by such corrections. Indeed, this perspective is consistent with an analysis 225j | that 


cons iders the effect of pion loops, which are explicitly excluded in the rainbow-ladder truncation 
H: whilst the individual masses are reduced by roughly 300 MeV, the mass difference, 5m, 
increases by only 50 MeV. Improving upon that analysis using insights drawn from the extensive 
body of work completed by the excited baryon analysis center (EBAC), which used a realistic 
coupled-channels mode l to remove meson dressing from the A and expose a dressed-quark-core 
bare-mass of 1.39 GeV mi, it is appropriate to choose 


gN = 1.18 , 5A = 1-56 ^ mw = 1.14 GeV, m-A = 1-39 GeV, 5m = 0.25 GeV . 


(3.7) 


This discussion introduces the basic elem ents in a contact-interaction computation of the 
spectrum of hadrons with strangeness 205|, which is the most extensive exploration of the 
baryon spectrum ever accomplished using the DSEs and yields numerous insights into hadron 
structure. Eor example, it predicts that the diquark content of baryons is largely independent 
of strangeness, viz. that baryon structure is flavour-blind; and the computed level ordering 
matches that of experiment, as apparent in Eig.[T2j In particular, the parity-partner for each 

tP - 1 - 
2 


ground-state is always more massive than its first radial excitation, i.e. the first ^ state 


always lies above the second st ate. Other approaches find this ordering difficult to 

achieve. That, too, was explained in Ref. |2n5| | : a veracious expression of DCSB in the meson 
spectrum is critical to obtaining the correct level ordering in the baryon sector, so that an 
approach within which DCSB cannot be realised or a simulation whose parameters are such 
that the importance of DCSB is suppressed will both necessarily have difficulty reproducing the 




















experimental ordering of levels. In addition, the contact interaction study in Ref. |2n5l | predicts 
that the first radial excitation of ground-state baryons is constituted almost entirely from axial- 
vector diquark correlations. A subsequent study [^, employing a QCD-linked interaction, 
predicts a completely different result; namely, that radial excitations possess the same diquark 
content as the ground-states. This is precisely the sort of thing one was looking for: signature 
empirical differences between predictions obtained with different interactions treated at the same 
level of approximation. 

Looking closer at Fig.[T2]it is evident that the computed baryon masses lie uniformly above 
the empirical values. This is a success because the computed masses are those of the baryons’ 
dressed-quark-cores, whereas the empirical values i nclude eff ects associated with meson-cloud 
effects, which typically produce sizable reductions 2261 . 227], as explained above. The values 
in Fig.[T2] may reasonably be viewed as bar e-ma ss inputs appropriate for dynamical coupled- 
channels analyses of the hadron spectrum 122811 . This was explained and illustrated for the 
nucleon and A-resonance in Sect. 4.5 of Ref. jlSa j and in particular for the Roper resonance in 
Ref. [HU. Here it is worth reiterating those instances in which a comparison can be made: 


(3.8) 


where when it appears, is the relevant bare mass inferred in the associated coupled- 

channels (DCC) analy sis. The DCC bare masses were uncertain and dependent on model details. 
However, as Ref. [205 ] made no attempt to fit them, their proximity to the DSE results suggests 
that it might now be possible to place these DCC bare masses on firmer ground, investing 
them with meaning within the context of hadron structure calculations that have a traceable 
connection with QCD. This is the subject of an ongoing investigation. 
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3.5. Predictions obtained with a Realistie Interaction 

The ultimate goal is ah initio predictions of hadron observables within continuum QCD; 
and for th e nucleon , A-baryon and Roper resonance, studies of the Faddeev equation exist 
57, 58, 64, I72I . I85I ] that are based on the one-loop renormalisation-group-improved interaction 
which was used efficaciously in the study of mesons. All these analyses retain the scalar and 
axial-vector diquark correlations that are known to be necessary and sufficient for a reliable 
description of positive-parity baryons. In order to compute baryon elastic form factors and the 
form factors describing transitions between them, it is necessary to know the man ner through 
which a composite nucleon couples to a photon. That is described in Ref. [229(1 ; and using 
this current one obtains predictions for all four nucleon elastic form factors, which are depicted 
in Fig.[T3j The notation is explained by recalling that the Poincare-covariant electromagnetic 
current for a spin-half nucleon is: 

J^{K,Q) = ieu{Pf) K^{K,Q)u{Pi) = ieu{Pf) , 

(3.9) 

where Pi (Pf) is the momentum of the incoming (outgoing) nucleon, Q = Pf — Pi is the 
photon momentum, K = (Pi + Pf)/2 is the total momentum of the system, K ■ Q = 0 and 
= —m|f(l -|- Tat), tn = (5^/(4m|^) for elastic scattering; and Fi and F 2 are, respectively, 
the Dirac and Pauli form factors, from which one obtains the nucleon’s electric and magnetic 
(Sachs) form factors 




F2(QQ, GMiQn = FiiQ^) + F2{Q^) . 


Ge(Q2) = Fi(Q2) 


(3.10) 


































Figure 13. Proton (top) and neutron (bottom) electromagnetic form factors. In both rows: left 
panel - Sachs electric; right panel - Sachs magnetic. Curves in all panels: solid, black - result 
obtained in Ref. [b^, using a Faddeev equation kernel and interaction vertices that possess QCD- 
like momentum dependence, viz. a QCD-kindred framework; dotted, blue - result obtained with 
a symmetry preserving treatment of a cont act in teraction (Cl framework) [^; dot-dashed,red - 
2004 parametrisation of experimental data 2301 ]. 


It is apparent in Fig.llSIthat the QCD-kindred results are in fair agreement with experiment, 
which is represented by the yea r-20 04 parametrisation in Ref. [23Cll |. (Comparisons made with 
a more recent parametrisation 23l| are not materially different.) No parameters were tuned 
in order to achieve this outcome. The most notable mismatch appears to be in the description 
of the neutron electric form factor at low Q^. However appearances are somewhat deceiving 
in this case because, on the low-Q^ domain: Gfj is small and hence slight differences appear 

such as 
bill, form 


large; and is much affected by subdominant effects that were neglected in Ref. 
so-called meson-cloud contributions. On the other hand, as was previously observed 
factors obtained via a symmetry-preserving DSE treatment of a contact-interaction are typically 
too hard. The defects of a contact-interaction are expressed with greatest force in the neutron 
electric form factor. One thus arrives at a completely unambiguous conclusion: 



Direct comparisons between experiment and sensibly formulated theory can distinguish 
between the momentum dependenee of the interaction that underlies strong-interaetion 
dynamics. 


In this connection one should guard against models with no traceable connection to QCD, which 
are tuned to fit data, and practitioners inclined to obscure this empirical fact. 

With predictions for the nucleon form factors in hand, one may revisit the empirical discovery 

















Figure 14. Left panel: normalised ratio of proton electric and magnetic form factors. Curves: 
solid, black - result obtained in Ref. 6^ using the QCD-kindred framework; Dashed, blue - 
Cl re sult (^ : and dot-dashed, red - ratio inferred from 2004 parametrisation of experimental 
data 230]. Data: blue circles |^ : green squares |^; brown triangles [s^; purple asterisk (^ : 
and orange diamonds (^ . Right panel: normalised ratio of ne utron electric and m agne tic form 
factors. Curves: same as in left panel. Data: blue circles mi; and green squares 


described in connection with Eq. (|2.2p . so in Fig.[T3]I depict the unit-normalised ratio of Sachs 
electric and magnetic form factors for the proton and neutron. Let us first consider the proton’s 
ratio (left panel). Both the contact-interaction and QCD-kindred frameworks predict a zero 
in but the comparison with extant experimental results indicates that the contact- 

interaction is invalid on > M(0)^, where M{p^) is the dressed-quark mass function explained 
in Sec. 12.51 The result obtained with the QCD-kindred framework, on the other hand, agrees 
with available data and predicts a zero in this ratio at ~ 9.5 GeV^. Notably, owing to the 
presence of strong diquark correlations, the singly-represented d-quark is usually sequestered 
inside a soft (scalar) diquark correlation. The appearance of a zero is therefore driven primarily 
by the contribution to G^ from the doubly-represented n-quark [^, which is four times more 
likely than the d-quark to be involved in a hard interaction. (Additional novel insights are 
provided in Ref. [65[.) 

As explained in Refs. j6l|,[6^, the behaviour of the dressed-quark contributions to the proton’s 
electric form factor on > 5 GeV^, and hence G^ itself, are particularly sensitive to the rate 
at which the dressed-quark mass runs from the nonperturbative into the perturbative domain 
of QCD. This is readily explicated using the information in the left panel of Fig. 1141 

The contact-interaction produces a momentum-independent dressed-quark mass; and in this 
counterpoint to QGD the dressed-quarks produce hard Dirac and Pauli form factors, which yield 
a rat io u ^Ge/Gm that possesses a zero at < 4 GeV^. Alternatively, the mass function used in 
Ref. [6^ is large at infrared momenta and approaches the current-quark mass as dressed-quark’s 
momentum increases. Such is the behaviour in QCD: dressed-quarks are massive in the infrared 
but become parton-like in the ultraviolet, characterised thereupon by a mass function that is 
modulated by the current-quark mass. Hence, the proton’s dressed-quarks possess constituent- 
quark-like masses at small momenta and thus have a large anomalous magnetic moment on this 
domain. As the momentum transfer grows, the structure of the integrands in the computation of 
the elastic form factors ensures that the dressed-quark mass functions are increasingly sampled 
within the domain upon which the transition from nonperturbative to perturbative behaviour 
takes place. This corresponds empirically to momentum transfers > 5GeV^. The rate 
at which the transition occurs determines how quickly the dressed-quarks become parton-like. 
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Figure 15. Left panel. Dressed-quark mass function employed in Ref. [6 
form, obtained in a least-squares fit to a diverse array of pion properties [23; 
the domain upon which DCSB is active. Right panel. Response of PpGe/Gm to increasing a; i.e., to an 
increasingly rapid transition between constituent- and parton-like behaviour of the dressed-quarks. Data 
are from Refs. [51, 53 - 5^ . 


i.e. how rapidly they are unclothed and come to behave as light-fe rmio n degrees of freedom. 
Since light-quarks must have a small anomalous magnetic moment 


then this transition 

entails that the proton Pauli form factor, generated dynamically therewith, drops to zero. This 
produces an interplay between the Dirac and Pauli form factors which, via Eq. (13.101) . entails 
that a momentum-dependent mass-function must beget a zero at larger values of than is 
obtained with momentum-independent dressed-quark masses. 

The dressed-quark mass function used in Ref. is characterised by a particular rate of 
transition between the nonperturbative and perturbative domains. If one were to increase that 
rate, then the transformation to partonlike quarks would become more rapid and hence the 
proton’s Pauli form factor would drop even more rapidly to zero. In this case the quark angular 
momentum correlations, expressed by the diquark structure, remain but the individual dressed- 
quark magnetic moments diminish markedly. Thus a more rapid transition pushes the zero in 
RpGep/Gmp to larger values of Q^. There is a rate of transformation beyond which the zero 


disappears completely [62l . l65l |: as made plain in Fig. 1151 there is no zero at all in a theory in 


which the mass-function rapidly becomes partonic. 

It follows that the possible existence and location of the zero in the ratio of proton elastic 
form factors [iipG^^{Q'^)/G^j^{Q‘^)] are a fairly direct measure of the nature of the quark- 
quark interaction in the Stan dard Model. Like the dilation of the meson valence-quark 
parton distribution amplitudes |l45l4l4^ . 235| |. they are a cumulative gauge of the momentum 
dependence of the interaction, the transition between the associated theory’s nonperturbative 
and perturbative domains, and the width of that domain. Hence, in extending experimental 
measurements of this ratio, and thereby the proton’s charge form factor, to larger momentum 
transfers, i.e. in reliably determining the proton’s charge distribution, there is an extraordinary 
opportunity for a constructive dialogue between experiment and theory. That feedback will 
assist greatly with contemporary efforts to reveal the character of the strongly interacting part 
of the Standard Model and its emergent phenomena. 

Let us return now to the right panel of Fig.[TH which displays the ratio The 

neutron ratio also exhibits a zero but at ~ 12GeV^, i.e. shifted to a 25% larger value of 
compared with the zero in the proton ratio. The properties of the dressed-quark propagators and 
bound-state amplitudes which influence the appearance of a zero in are qualitatively 

the same as those described in connection with fipG^/G^. However, owing to the different 
electric charge weightings of the quark contributions in the neutron, the quantitative effect is 



























Figure 16. Left panel. Ratio pnG^/G% for N = neutron, proton: solid, black - neutron result obtained 
in Ref. [^ . using QCD-like momentum-dependent quark dressing; dashed, blue - neutron result obtained 
with such dressing but an accelerated rate of transition from dressed-quark to parton, viz. a = 2; and dot- 
dashed, red - proton result obtained in Ref. (^ . dotted-green - proton result with accelerated transition 
{a = 2). (The arrows highlight the response to accelerating the dressed-quark—>■ parton transformation.) 
Right panel - Ratio of neutron and proton Sachs electric form factors: solid-black - result obtai ned in 
Ref. [bj; and dot-dashed-red - ratio inferred from 2004 parametrisation of experimental data |230l| . 


opposite to that for the proton. Namely, when the transformation from dressed-quark to parton 
is accelerated, as described in Ref. and illustrated in Fig. 1151 the zero occurs at smaller Q^. 
This is depicted in the left panel of Fig.1161 On the other hand, as indicated by the dotted curve 
in the lower-left panel of Fig . 1131 one typically finds that a contact interaction produces no zero 
in the neutron ratio 


61 


The origin of these features is readily elucid ated . Note fir st that the s-quark contributes very 
little to nucleon electromagnetic form factors 223l . 2361 - 238 ] and therefore write 


— BuG^^ — \eii\G^'^ , 

where the isolated terms denote the contribution from each quark flavour, 
charge symmetry is almost exact in QCD, so that 


(3.11) 

Consider next that 


^n,d _ ^P,TJ 
'^E ~ '^E 




_ ^p,d 

— '^E ■ 


(3.12) 


and hence, to a very good level of approximation, 

Gg = BuG^^ — je^lG^’^ = BuG^^ — |ed|G^“ . (3.13) 

Now, with a zero in G^ at k. 9.5 GeV^ =: Sz, one has G^'^{sz) = 2G^“(s^) and hence 
G%{sz) = G^^{sz) > 0, where the last result is evident in Fig. 7.3 of Ref. [^, which shows that 
although the behaviour of G^“ and G^ is qualitatively similar, the zero in Gg“ occurs at a 
larger value of than that in G^ itself. Under these conditions, any zero in Gg must occur at 
a larger value of than the zero in G^: compare the dot-dashed and solid curves in Fig.fTBl 
This relative ordering of zeros can change, however, because, in contrast to Gg“, G^^^ evolves 
more slowly with changes in the rate at which the dressed-quark mass function transits from 
the nonperturbative to the perturbative domain, something which is also apparent in Fig. 7.3 of 
Ref. [^. As explained above, this inertia owes to the d-quark being preferentially sequestered 
inside a soft (scalar) diquark correlation. Subject to these insights, consider Eq. (13.13^ : with the 











Figu re 17. Comparison between CLAS data 
|239l | on the magnetic 7 + TV —>■ A transition 
form factor and a recent theoretical prediction (^ . 
The dashed curve shows the result that would be 
obtained if the interaction between quarks in QCD 
were momentum-independent The solid curve 

is obtained with precisely the same QCD-based 
formulation as was employed in a successful analysis 
of nucleon elastic form factors, which explains 
the behaviour of the ratio /rpG^(Q^)/G'^(Q^) 
described in connection with Eq. (12.21) and elucidated 
above. The experiment-theory comparison confirms 
that experiments are sensitive to the momentum 
dependence of the running couplings and masses in 
QCD; and the theoretical unification of N and A 
properties highlights the material progress that has 
been made in constraining the long-range behaviour 
of these fundamental quantities. (Figure courtesy of 
V. Mokeev.) 


location of a zero in shifting slowly to larger values of but that in moving rapidly, 
one is subtracting from a function whose domain of positive support is becoming 

increasingly large. That operation will typically shift the zero in to smaller values of and 
eventually enable a zero in even when that in the has disappeared. 

The right panel in Fig.[T 6 ]displays a curious effect arising from the faster-than-dipole decrease 
of the proton’s electric form factor (and possible appearance of a zero); namely, there will likely 
be a domain of upon which the magnitude of the neutron’s electric form factor exceeds that 
of the proton’s. This being the case, then at some value of momentum transfer the electric 
form factor of the neutral composite fermion becomes larger than that of its positively charged 
counterpart. That occurs at = 4.8M^ in the QCD-kindred analysis of Ref. [h^. JLab 12 will 
test this prediction. 

I would like to close this section by emphasising that given the challenges posed by non- 
perturbative QCD, it is insufficient to study hadron ground-states and elastic form factors alone. 
Many novel perspectives and additional insights are provided by nucleon-to-resonance transition 
form factors, whose behaviour at large momentum transfers ca n re veal much about the long- 
range behaviour of the interactions between quarks and gluons 240l |. Indeed, the properties of 
nucleon resonances are more sensitive to long-range effects in QCD than are those of hadron 
ground states. The lightest baryon resonances are the A(1232)-states; and despite possessing 
a width of 120 MeV, these states are well isolated from other nucleon excitations. Hence the 
'y + N ^ A transition form factors have long been used to probe strong interaction dynamics. 
They excite keen interest because of their use in probing, inter alia, the role that reson ance 
electroproduction experiments can play in exposing non-perturbative features of QCD 2401]: th e 
relevance of perturbative QCD to p roce sses involving moderate momentum transfers 241142431 ]: 
and shape deformation of hadrons |244l |. Using the “CLAS” detector at JLab, precise data on 
the dominant 7 -|- A —>• A magnetic transition now reaches to = 8 GeV^; an eventuality that 
poses both great opportunities and challenges for QCD theory, some of which have recently been 
met, as illustrated in Fig.fTTl 













4. Epilogue 

This series of lectures has focused on quantum chromodynamics (QCD), the most interesting 
part of the Standard Model and Nature’s only example of an essentially nonperturbative 
fundamental theory. They have introduced the Dyson-Schwinger equations (DSEs), which 
provide a nonperturbative tool that may be used to gain insights and work toward a solution of 
QCD in the continuum. In progressing through these lectures we have seen that the DSEs have 
enabled the proof of numerous exact results in QCD, amongst them: 

• quarks are not simply quantum mechanical Dirac particles, they are far more complex 
objects; 

• gluons are nonperturbatively massive; 

• dynamical chiral symmetry breaking (DCSB) is a fact - its responsible for 98% of the mass 
of visible matter in the Universe; 

• Goldstone’s theorem is fundamentally an expression of near-equivalence between the one- 
body problem and the two-body problem in the pseudoscalar meson channel; 

• confinement in our real-world is a dynamical phenomenon, which can neither be expressed 
via any standard quantum mechanical potential nor discussed sensibly in the absence of 
dynamical light-quark degrees-of-freedom; 

• condensates, those quantities that have commonly been viewed historically as constant 
empirical mass-scales that fill all spacetime, are instead wholly contained within hadrons, 
ie. they are a property of hadrons themselves and expressed, for example, in their Bethe- 
Salpeter or light-front wave functions; 

• etc. 

Moreover, we have seen the application of DSEs in explaining and predicting a wide range of 
experimental data. Thus, at the conclusion of these lectures it will be clear that the DSEs are a 
single framework which is almost unique in providing an unambiguous path from a defined 
interaction —)• confinement and DCSB —)■ hadron masses, radii, elastic and transition form 
factors, distribution functions, etc. Namely, the DSEs are a practical, predictive, unifying tool 
for fundamental physics. 
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